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The following pages describe the solutions to many of the Exercises in 
David Chandler's text. Introduction to Modern Statistical Mechanics . That 
text is brief but demanding with the principles illustrated only through the 
numerous Exercises. Especially in light of this strategy, a solution manual 
can serve an important pedagogical role. While we do not present herein the 

1P>' 

solutions to all the Exercises, those we have skipped require no techniques 
that are not already illustrated by the solutions we have included. In some 
cases, we have pointed out additional references. While not absolutely 
necessary (the text is self contained), these additions may be especially 
useful for the topics not covered explicitly in the main body of the text. 

For instance, in Chapter 6 where at least one exercise leads one naturally to 
perform a finite size scaling calculation, we alert the reader to some 
introductory literature on this important technique. 

In the format we have adopted, the manual should look like neatly 
prepared problem sets. The prose is kept to a minimum, and the figures are 
generally hand drawn sketches, though a few graphs were plotted with the aid 
of a microcomputer. Concerning the preparation of the manuscript, we owe 
great thanks to Mary Hammond for her expert secretarial help. 
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Chapter 1. Thermodynamics, Fundamentals 



1.2 Entropy is postulated by the 2nd Law to be a monotonically increasing 
function of E (i.e., T > 0) and an extensive function. The first 
equation of state satisfies both criteria, while the second is not 
extensive since it grows exponentially with system size. So 



S - L 0 Y(8E/L 0 ) 1/2 - L q Y[ 1 [k-f ♦ ^ - | ] 



is the correct choice. Indeed, you can convince yourself that 



S( AE, AL, An) = AS(E,L,n) . 



Further , 



1 



1 - fiih = 1 
T 3E n,L 2 



implying T > 0 . For the tension, we have 



-f/j - = - ly [ — - — 1 

1 c) L E , n L 0 L 2 . 2 J 

L 0 L 
3 

or f = TY( Sl/dg ) [ 1 - (Hq/S.) 1 , Z = L/n . 



1 . 4 The principles we apply are conservation of total extensive energy and 
equality of temperature at thermal equilibrium. From the solution to 
Exercise 1.2, 

e = t 2 y 2 ni.e/'i . 

■Rf&: v 



Therefore, the total initial energy is 



[(T (1) ) 2 n (1) + (T (2) ) 2 n (2) ] Y 2 R. o 0/-4 = E (1) + E (2) 
which is the same as the total final energy 



T 2 (n (1) + n (2) ) Y 2 a o 0/i4 . 



1 



T - ^[(r (,) ) 2 n^ + (t< 2 >)V 2 > 3 l 1 ' 2 



ijg 



SdT + Ldf + ndp 



1.13 Since, 



1.5 Assuming the system is surrounded by adiabatic walls, the change in 

energy AH is equal to the work jpdV due to changing the pressures p ft and 
p B . The final equilibrium state so obtained has new mole numbers n, 
and n 2 • Instead we can imagine arriving at this final equilibrium 
state in two steps. First we reversibly pump the particles to their 
final values. Then we do work by moving the pistons to their final 
pressures with the walls impermeable to any particles. We know the 
total change in energy of the system must again be AE since E is a 
function of state. Hence, the work done by changing n, and n 2 can be 
associated with work done by changing the pressures. 



1.1 



T (1- (!§.) 

1 'JT'r 



3C 

^3p~^T,n - 1 '■Sp '3T'p,n'T,n 

. T f— ( M) ) 

*-31 3p T,n J p,n 



T ri_ ( - M) } 

U 3T 3T p,n J p,n 

- T(— ] 

T 3T 2 P ,n 



p,T,n natural variables => 
dG - -SdT + Vdp + pdn 



! 

i 

I 

§ 

j 



2 2 
E = eS L/n 



we have 



f 3E l 

u 3n S, L 



20S 2 L 



Substitute S = S(T,L,n) by inverting T = T(S,L,n): 
t = f 3 !] 2eSL 



" 3S ; L,N 



n 



2 ’ S0 S = M ' 



Thus , 



y = 



29 ^) 2 1* 



n 



20(L/n) - 



The Gibbs-Duhem like equation is 0 = SdT + Ldf + ndp . We want to 
show this interdependence of these three differentials. We choose a 
representation, e.g., |s,L,n}, though we could also use { T , f } : 



dT 



^3S^L,n dS + ^3L^S,n dL + ^sJs.L dn 



= (^)dS + (20|) dL + (- 



40SL- 



dn 



1.12 For a rubber band of length L, tension f, 
dE = TdS + fdL + pdn . 

E is extensive. Therefore by Euler's theorem, E = TS + fL + pn 
Hence , 



Similarly for df and dp 

f . [Mi - 0 sf 

3L S,n 0 2 

n 



df = (£§§)«* + (- 



20S 



)dn 



du = (- ^)dS + (- — )d L + (i®|^)dn 
n 3 n 3 n 



dE = TdS + SdT + fdL + Ldf + pdn + ndp 
implying 



2 



3 



n 



Plugging into 0 = SdT + Ldf + ndp 



[SI— J + L[— + n[- 5— J ] dS 

n n n 3 

+ CS(-^|) + L(0) + n(- ^|-)] dL 
n rr 

r „f 46SL , ( 20S 2 i fbes 2 !,-, . 

+ [S[- — y- + l[- — 5 -J + n(— — J dn 

n J n- 3 n 



1 . 1 M The Gibbs-Duhem equation, 0 = -SdT + Vdp - ndy, implies 
dy = -sdT + vdp , with s = S/n and v = V/n . 



Hence , 



(— ) = 
'-Sv-'T 



^T + V ^T 



vf — ) . 

'■3v'T 



1.1 6 From the definition of heat capacity 



f_i) . ) . t( — f— ) ) . 

1 3J> J T l 3S, U 3T J 8/T *■ 3T 1 - 3il J T'Jl 



Further, since dE = TdS + fdl + ydn, we have 

d(E-TS) = - SdT + fdL + ydn 



and thus the Maxwell-type relation 



r 3s -I r 3f i 

L 3fc J T = " 3T H 



3o. a 2 f 

[— ] = - Tf— ] 

d I 



t 4 (?) ■ ° 

3T 






are 



Chapter 2. Conditions for Equilibrium and Stability 



01 . 

(a) Stability criteria are generally (^i) > 0 where X, and I- 

0 A “* 11 

1 



conjugate. Therefore (i) is true, 



= n ^fv^T n ^ ^ ’ since “P an d V are conjugate. 

Also (iii ) is true since (|H) t = + v (f^ T < 0 by the Gibbs-Duhem 
equation. None of the remainder are guaranteed false by stability. 



(b) From Part (a), we know to examine (ii) and (iv). Note 



fil) = - fill flil 

l 3v J s l 3s'v l '3v J T 



(— ) (—) 



' 3s ; v l 8T J v ‘ 



{dA = -SdT - pdv + ydN} 



Since (3T/3s) v is positive by stability, 



> 0 implies (|£] < 0 . 



But this implication contradicts (ii). Thus, (ii) and (iv) are 
inconsistent . 



2.10 To derive the analog of the Clausius-Clapeyron equation in the y-T 
plane, note that the y-T coexistence line requires that 
p^ a hy,T) = p^(y,T) . But 0 = -SdT + Vdp - ndy implies 
dp + (S/V )dT + (n/V)dy . As we creep infinitesimally along the y-T 

(a) ( a ) 

coexistence line, dp = dp when infinitesimally close on either 
side, i.e., 



(|) (a) dT + (tt) Cct) dy - (|) (B) dT + (B)(6>d u 



|-n-, (a)_ (- n -> (g) 

Ap - l v J lyj 

|-S-j(a) (- S> ( B ) 

As = - L^J 



2-!9 C v -T(|f) v 



c - t(— ] 

p 1 3T p 



1 (31) 

V '■3p J s 



K - - IfMl 

T V'-Sp-'T 



Assume constant n . To prove (K g /K T ) = (C y /C ) we note that 

V k t = 

f 9V i r9S-i f3pi r3T'| 

'■3S J p l '3p' l v l '3T J v''3V J p 



_ (M] r 3Tl _ (3§1 /(M) = JL _ 

l 3T J v l '3S' , p '■3T J v '■3T J p C p 

, . . ( 3xi ( 3xi r 3 z -v 

The second equality comes from realizing lgyJ z = ~ '■3z-'y''3y-'x 
and the third uses the chain rule. Since 



Thus , 



c - c = - t (£E) r (_) ^ 

p v 1 3V T,n L 3T p ,n J 

and since -Qp/3v) p > 0 by stability, we have C p > C y 
C p /C v >1, and since C p /C v = K ? /K s , we have K ? > K g . 

This inequality says that instead of working adiabatically , it is easier 
to compress something if you can leak energy out of the system into a 
constant temperature heat bath. 



2.20 (a) The rubber band heats up when stretched. Assume constant n. So, 

©s,^ 0 0r e ^ivalently © s>n > 0 . 

Note, these two derivatives have the same sign since 



r 3T 



(— ) f— ) 



3L i 



'■SrS.n '-3L ; S,n v 3f'S,n 

and the second derivative is positive by stability. 

To find the sign of (M) , we write 

ol i >n 



2.21(a) 



mm 



r 3L ^ ^ 9L 'j ( 9S -j 

^ r , n " l 9S j f,n 



'■3f ; S '■3T J f 



f 3Si 



>0 >0 

given by stability 



d(E-fL) = TdS - Ldf + ... 

• ,. f 3L -> r 3T i 

implies (gg) fjn = - ( w ) s>n 



Thus 



'■3T J f ,n 



< 0 . 

In other words, the rubber band stretches when cooled. 



(b) 



(M) = fi£) fii) + (M) (9b) 

(. 9T J f '■axJL l '3T J f ^3L J T '■3T J f 



- r^] = f^] r-f— 1 f— 1 i 

3T ; f l '3Ta '■3L- , T l ^ 3T'L v 3f'T J 



r3Si 



r 3S-, 



r 9f 1 



r 3L 



= (M) (M) 

‘-OL-'t ^ 3f -'t 



' 3L-i 



d(E - TS) 

- -SdT + f dL + . . . 



>0 



> 0 

by stability 



implies 

(H) , 

l '3T J L 



l 3L J T 



c L > 0 



Hence 

L 'L f 



In other words, the constant length rubber band 
has the larger change in temperature. 



Since in general dE = TdS + f*dx , where f is the intensive field 
and x is the extensive parameter, the correct work term is HdM. 
(In other words, doubling the system doubles M, not H. This 
approximate extensivity neglects interactions such as mutual 
polarization which might occur if two magnetizable systems were 
placed adjacent to each other.) So 

dE = TdS + HdM - pdV + pdN , 
and 

dA = -SdT + HdM - pdV - pdN . 



Stability then requires both x 1 = (Mb ,, > 0 and 



k 3M'S,V,n 



-1 _ r 3H-i 

X T l 3M J T,V,n 

Using the representation {s,H,n,v} , the thermodynamic potential 
E - HM. We then have M = M(S,H,n,V) since M = { — ^S.n.V ' 

f 3M i _ r 3M i f3Sl 
X T " l 8H J T ^ 3S H ' ' 3H T ^ 3H S 



C 3T-) t STi . r 3S-) 

X T ' X S r 3H J S *- 3H S '■3T J H 



= fill /fill >0 

l 3H ; S l '3S J H 



d(E-HM) = TdS - MdH + .. 
By the Maxwell relation 

rMi = -fill . 

' 3S H V 3H J S 



(||) H > 0 by stability. 



(b) Given (|=) H < 0 , we want to determine the sign of (|^) s = 



fill = - fill .fiil 

'■3H j S ''SS'H ^SH^T 



d( E-HM- ST) = -SdT - MdH + 



= - fill fiih „ 

^3S"n l 3T J H 

"■ — ? — 1 ' — t — ' 

>0, by <0, 
stability given 



By the Maxwell relation 

flil . (Ml . 

1 -3H ; T ^3T ; H 



2.23 (a) dA = -SdT + fdL + pdM 



x = L/M 



k, h, x n and c are independent of x, but not T. 



f S " [ 3 eHt,M - MkX- M ■ kX 



f L = Mh(x-x Q )- R = h(x-x Q ) 



(b) A - 1 kMx 2 = 1 kL 2 /M 



u - (— ) 

U S k 3M 'T,L 



L 1,2 

— = - - kx 

,2 2 



A l - \ hM(g - x Q ) 2 ♦ CM - 1 h[£- - 2Lx q + Mx 2 ] + CM 



Since x n is independent of L and of x, it's independent of M. 



p L = i h[ - 7 + x 0 ] + 0 - - i hCx2 - x o } + c 

M 

(c) - f s x = - 1 kx 2 = Ug 



L f „ 1 . , 2 2. 

r - f L x = - - h ( x - x 0 ) + c = y L 



(d) The phase transition occurs at u = u - n f - r r 

M S U L U T , Ig - t L - f T . 

At p T> f T , the short phase has x = x g , and the long phase has 



P s = p L implies - | x | = - | (x 2 - x 2 ) + c . 



f S = f L iraplies kx c = h ( x, — x„) . 



S L 0' 



(kXg) 2 = hkx 2 - hkx 2 - 2ck , 



X L - X 0 ” 



hkx Q ± /hk Akx Q + 2c(k-h) 
h(h-k) ~ 



f T h(x L X 0 ) 



since k>h, c>0 , we have, 



/hk /hkx^+ 2c ( k-h) - hkx Q 



^ x s - £ (X L - x o^ iraplies Ax = x - x 



HH ( V V + x o 



Ax - / h/k /hkx Q + 2c (k-h) - hx^ + x^ 



•2iJ Near the liquid-solid phase coexistence 



, (£) 1 a n c 



2 T V 



,(s) 1 1 

3 T 2 
j V 



which are functions in (n,V,T} . 



( fif ) ( s ) 

(a) p and p are the densities near the coexistence line. At 
coexistence, the temperature is the remaining free parameter, 
P (S) - p U) , p (3) = p a) , T (S) = T U) - T . 



Explicitly, 

a , (a) a .(a> 

(M. ) = 1 

'■an J V,T '■an J V,T 



'3A (s \ 

' 3V J n,T ’■3V J n,T 



2 ( 0 ) 

irf-) 

T l ' v J Tv' 



2 B r /flv (s) = a r ,ru (ft) 

3 X L ''v ; J 2 T LV v j 



8 (s) * (5,) 

Substituting — [p ] = p into the previous equation we get 

r (s) 3 3a r g (s) 2 2 
Cp ] = H6 C a p ] ’ 



Therefore the densities are 



(s) b a . (Jl) 16 a 
p =33 and p - ^ - , 



and are constant functions of T. 



(c) S k 



( i) a 


r < «•) >2 


- P 2T 


Ip J = 


3A (8,) i 


2 

1 a n 


■ 3T n,V 


2 t 2 V 


As) _U) 

: “ 0 


1 r 0 


n 


T 2 '3 P 



81 „2 T 



(s)_ - + 1 i_ HI 

’ _ l 3T J n,V 3 t 2 v 2 



1 ) - i_r_ — - - 
1 " t 2 L 27 6 9 8 



where s = S/n 



(d) Since Av = 



1 6 a ’ 



As = J_28 or 1_ 

AV " " 81 8 2 T 2 ' 



dp s = _ _1_28 a 3 1_ 
dT 81 g 2 t 2 ' 



Therefore _ _ , satisfying the Clausius-Clapeyron equation. 



.25 Stability requires - (|£) > 0 

dv i ,n 

or flEl - - if flEl > 0 
3p T,n n l 3V J T,n > 0 



From the van der Waals equation, 



f9P] _ RT 
3p T,n (1-bp) 2 



2 ap 



Instability occurs when this function is negative. The boundary 
surrounding the region of instability, the spinodal, is where it’s 0: 



RT - 2ap( 1 -bp J = 0 . 



In the T-p plane, the spinodal is RT = 2ap(1-bp) 2 , i.e., the spinodal 



is a cubic polynomial with roots for T = 0 at p = 0, 1/b . Also note 
that the spinodal has extrema when 



p p 

= 0 = 3b p - bbp + 1 , 



that is when p = 1/b , 1/3b . 



Similarly, there is a point of inflection when 
2 

d T 3 

— - = 0 = 3b J p - 2b = 0 , 

8p 

that is when p = 2/3b . 



Finally, 



Tio 1 n _ 8 a 

p 3b 27 b 



These features are illustrated in the figure below. Only p < b 1 is 



considered since higher p would give negative pressures. Note also that 



the pressure on the spinodal 



Van der Waals Isotherms and Spinodal 




becomes 



- PtT anlnft , nl (p).p3 - a e (1 ' 2bp) 



^spinodal H spinodal 



negative for p > 1 /2b 



Van der Waals Spinodal 



2.26 phase a: gp = a + bgp 



phase g: gp = e + d(gp)‘ 



At phase equilibrium, 



The van der Waals equation of state for a given temperature T Q has 
only one interval in p, if at all. for which (|£) T < 0 , i.e., for which 
the fluid is unstable. This is reflected by the fact that the spinodal 
equation, p = ap 2 (1-2bp) , has only two distinct roots. 

in other words, the isotherm is shaped as illustrated in the figure 
on the next page. The Maxwell construction requires p(p g ) - P<P t > which 
can only be satisfied within the height h. Any pair of points so picked 

enveloping the spinodal region. 



Thus a + bgy = c + d (Bu) , implying 



The density can be obtained from the Gibbs-Duhem equation, 

, „ , \ ( ( 3 ( Bp) i 

dy = -sdT + vdp , i.e., (1/v) = p - ^ = gJjJp * 



will lead to a coexistence curve 




and vie identify the positive root as the physical root to the quadratic 
equation. Hence 




Chapter 3. Statistical Mechanics 



The probability in a canonical ensemble is 
P(E) « £2(E)e' eE = exp[StnS5(E)-SE] . 

We assume that P(E) is so sharply peaked that <E> occurs at the top (or 
very close to the top) of the peak. The following steepest descent 
calculation will validate our assumption: Expand 8,nP(E) = Hnq(E) - BE 

about its maximum value 2,nP(<E>) which occurs at <E> , i.e., 



S.nP(E) = 8,nP( <E> ) + 6E 38, ^ (E) | <E> - j 



1 2 3_tnP(E) , 4 

2 I <e> 



0, because S,nP(E) 
has a maximum at <E> 



BE 



negative 



B 2 StnP, 



3 (Enfl-BE), 



3E 



2 l<E> 



3 S,nP j 

1 CE> = 3E 2 '<E> 






3 j- 3S,nQ( <E> ) ' 



3E 

r 3 . 3&nq 

1 '3E l_ 3E _ '-'<E> _ 3<E> 1 3<E> 

1 3T 1 

" k t 2 3<E> = “ k t 2 c 
k B T B 1 U V 



3B 

3<E> 



Hence 



P(E) a exp[- (E ~ < § > — ] 

2 V S 

For 6E = 10 -6 <E> , N = 1 0 21 , using <E> = 3/2 Nk g T , = 3/2 Nk 0 , 

-1 o 9 

P(E)/P(<E>) = e » 0.000 0001 



a billion zeros 



That's improbable! 



9 m/N = 



1 + e 



Be 



implies 



1 



N _ 

N-m ” 1-m/N 



= e - ge + 1 , 
„Be 



and 

N-m 

m 



= N/m - 1 = e 



Be 




So 



N! 



S/k = S,nn(E,N) = inf-rr; r-; r 

B (N-m)!m! • 



where ra = m(E) = E/e 



= N S,n N - (N-m) Zn (N-m) - ra Sin m 



„ „ r N i „ rN-nu 

= N Sin (- — J + ra ta [ J 

N-nr v m 1 



N In (e' Be + 1 ) 



Clearly, Slim S/k = 0 . 



1 + e 



g£ Be - S( B , N)/k g 



Also, 



S/k B (E.N) - N to {^j) ♦ m(E) tn - 1) 

= - N Sin (1 - 



3 (S/k 



So 1 /k T = 

D 



B' 



5—) + E/e Sin (g— _ 1 ) • 

1 



3E >N 



Sin (^| - 1 ) 
e E 



The last result is in agreement with the earlier formula, 

Qp 

N/m - 1 = e . Notice, however, that 1 /T can be negative when 

Ne , . , . E s e . , . 1 , - 

1 < 1 , i.e., jj > 2 implies - C ° . 



1 



1 




For 



E and N varying, the grand canonical ensemble yields 



.2 r 3<E> -i . 2 r 3<E> i 

<(6E) > - l 3 g J v ,By ~ k g T l 3T J Vigy 



i.; e could calculate this last derivative precisely (see below), but for 
the purposes of an estimate, we only need to know that we'd only get 
some constant times C v , the exact constant being irrelevant since it's 
still 0(N). So, as an estimate 

<(8E) 2 > = k B T 2 C v 

which is the canonical ensemble result. The RMS deviation 



^<(7e? 

<E> 



/k 0 T 2 3k B N/2 

3k D TN/2 

D 



1 -1 1 
— *10 

•/n 



Similarly, for an ideal gas, <(6N) > = <N> 

,-11 



SO 



/< 



'<(6p) 2 > 






6N) > 



<P> 



<N> 



1 

/N 



1 0 



s-3uat for reference, let's evaluate the actual derivatives involved in 
<(6E) 2 > : 



<■ 3 E ’j _ /■ ^_E -j (—1 

l 3rV,By ~ l ‘3T J V,n '■3N j V,T '■3T'V,By 



3.10 From p. 68, the canonical ensemble calculation gives 



-BA = N in(1 + e~ Be ) 



<E> = Ne( 1 + e g£ ) 1 



Then 

S/k 



B 



-B(A - <E>) = N Hn(1 + e~ Be ) + 



NBe 



.Be 



1 + e h 

which is the same as in Exercise 3-9 from the microcananical 



“ c v + (3/2 k B T) • (- If^V.N 7 ^V,T • 



flini _ 5 (Mi + - 

3T ; V,N N l '3T J V,N N 



i 



- E . (- . 

N ^ 2 

k B T 



_3 

2T 



S = 1/T dE + p/T dV - y/T dN 
0 = EdB + Vd(Bp) - Nd(Bu) 

=> d(By) - E/N dB + V/N d(Bp) 



calculation . 



and 



22 

3.16 N = 0.01 moles is = 10 . An open, thermally equilibrated system is 

controlled by T, V and y. Fluctuations in the energy and the density 
p = N/v are given by 



<(6x) > 



3<x> 

36 



where g is conjugate to x 



I - 9 Bu -j 
3N J V,T 



V 3(Bp) | 

N 3N ' V, T 



V J_ _ 1_ 
N ‘ V ” N 



BP = ^ = Bp(N, V , T) 



rdE-j 3 

oT V, gy " C v + ( 3/2 KgT)-^)^ - (3/2 - 9/4)Nk 0 = 15/4 Nk 0 . 






above theorem on 



17 a < x < b . Show LVg,f|<gf> = <gXf>] <=> P( x ) = • 

(a) 

b b r b 

<gf> - <g><f> = / g(x)f(x)p(x)dx - j g(x )p(x)dx • J f (y)p(y )dy 



b b r b , b 

= f J g(x)f(y)p(x)< 5 (y-x)dxdy - J J g(x)f(y )p(x)p(y)dxdy . 

a a 







E 

n 




MM 




if 




(a) 









b b 

n 

a a 

Hence, if <gf> = <g><f>, 
b b 

/ / g(x)f (y)p(x)[s(y-x) - 
a a 

This holds for all g and f 

p(x)[6(y-x) - p(y)] - o . 

The trivial solution, p(x) = 0, is discounted because it is not 
normalized. Thus, we must have 



p(y) = 6(y-x Q ) , 



where x Q is some point in the interval [a,b] . 

(b) <f(x)g(y)> - <f (x)Xg(y )> 

= j j dxdy f (x )g(y )p(x ,y ) - J dx f(x)p^(x) / dy g(y)P 2 (y) • 

Hence, if <fg> = <fXg> , 

/ J dxdy f (x)g(y)[p(x,y) - P^x) P 2 (y)] = 0 • 

This holds for all f and g if and only if 
p(x,y) = p 1 (x) p 2 (y) . 

Note: The above proofs rely upon the fact that jf(x)n(x)dx -Of 

all f is equivalent to n(x) = 0 on the interval of integration. j 
This is true because if n(x) * 0 we could pick f(x) = n(x) ior 
which ff(x)n(x)dx > 0 . For two dimensions, we can iterate the 



/dx f (x)[/dy g(y ) n(x,y)] = 0 
calling the term in the square brackets n(x) 



18 For N spins in a field H, 
N 



= 1 (-n pH) , n = ± 1 

r 2’ •• N i=1 



-BE 



Q - l 



n 1 ,n 2 ,...n N 



-B I (-n pH) 
i-1 



n l’ n 2*--- n N 



I 



n l ,n 2’ ' ' ' n N 



n 1 ,n 2 ,...n N i-1 



N Bn . pH N 

n e 1 = n 1 e 

i=1 n=±1 



BpHn 



since the factor for each spin is identical, so 
Q _ (e ^H + e ‘ 6,jH ) N = (2 cosh(BpH)) N . 



as 



<E> = 



3(>nQ 

3(-B) 

- NpH 



g -BpH _ e BpH 
e -BpH + e BpH 



NpH tanh(-BpH) 






' 

' 

X. 



- NpH tanh (BpH) 



(b) inQ ■ k- - ITT “ - 6A 

B B 



s = 



-A + E 



= kgilnQ + k B B<E> 

= Nk D Un(e SyH + e~ 6>jH ) - BpH tanh (BpH)] 

B 

(c) In the limit T -*■ 0 or 8 ■* » , tanh (BpH) -» 1 . So 



<E> 



T->0 



= - NpH 



i.e., the ground state has all spins + , or aligned with the field 
lira Nk D [S.n(e 6uH + e~ 8yH ) - 6 pH tanh(BpH)] - Nk [BpH - BpH] = 0 , 

„ „ b D ' V 

hence , 



Vo ■ 0 ■ 



N 

l 

i=1 

Np tanh(BpH) 



N i N 

3.19 (a) <M> = < l pn . > = -•£<! - pn.H> = - ^ 

1 H 1=1 1 H 



2 

(b) <(6M) 2 > = <(M-<M>) 2 > = <M 2 > - <M> 2 = 1 3 Q 



J_ ( 3Q )2 

Q „2 ^a(SH) J 



3(BH) 



1 3Q 72 ■ SarlntNpQtanhCBpH)] 



3 (BH) 



= Np 2 (('N-1 )tanh 2 (BpH) + 1 } 



Hence , 



< ( 6M) 2 > = Np 2 { 1 - tanh 2 (gpH) ] = Npsech 2 ( BpH) . 



Also , 

3<M> | 

3H 1 B ,N 

„ r3<M> 



= BNp 2 { 1 - tanh 2 ( BpH) } . 



w>» ■ ■ «« 2 > • 



(c) As B ■+ “ , tanh(BpH) -> 1 



Therefore <M>„ = Np , and <(<5M) 2 >„ = Np 2 (1-1) = 0 . 

| CD 1 CD 

In other words, the ground state with all the spins aligned has no 
fluctuations . 



3.20 In an ensemble with constant M, all members have the same number of 



spins up. For a given H, this ensemble is the microcanonical ensemble 
since E = - HM is also fixed. Note we can't generalize this ensemble to 
one with fluctuating extensive variables because there is only one: M, 

which we fix (assuming N is constant), i.e., we can't have E flowing 
independently of M from the system to an energy bath. 

So we want to use E, H, and N as our natural variables. 



In this ensemble, 



M = u I n 1 = p(n + - n_) 
i=1 



= p(2n -N) , 



where n + = number of spins up 

n_ = number of spins down 

and n + n = N 



and 



E = - Hp(2n -N) . 



Since 



N! 



0(E) = n ! (N-n )! and S/k B = Ann(E) ' 



B = 



So 



J_ IS 

k B 3E 



i 3inn(E) 3inP 3n + 8 f „ r . 

'n,H = 3E = ^iT 3F = - H l Kn[n + (N-n + ) ! ] ) • 



-2Hp 



8H “ % 3 R- *n[n + !(N-n + )!] . 



For N >> 0, 



BH = Ip 3n~ V* 10 n + + (N-n + Hn(N-n + ) - N] 



- [in n + - ln(N-n + )] = ^ in (jp±-) . 



Thus, 



+ 1 ♦ e" 28 ^ ’ 

hence 



H = p(2n + - N) = Np tanh(BpH) 



which agrees with Exercise 3-19 and (after a change in notation) with 



BA -BA 
= e + e 



3.9 as well. 

Substituting n = M/2 y + M/2 into the above equation for BH gives 



in the 



{ A , B } basis H 



0 1 
1 0 



»» ■ k l " ' 



1 



and e' BH - e 



BA 



0 1 
1 0 



- I 

n=0 



[BA 



0 1 
1 0 



n! 





0 


> 


and m = 


u 


o] 


. 




-A 


0 


0 








V 


J 




V 


j 


■ .. 



3.21 In the { A, B } basis, H r 



In Pauli spin matrices notation, H Q - - Aa x and m = po 2 



- I 

n=0 

even 



( BA) 



1 0 
0 1 



I 

n=0 

odd 



(BA) 

n! 



0 1 

1 0 
t j 



Then H = H. - mE = -Act - pEo 
u x z 



(a) E = 0 



H = H 0 = - Aa x = 



0 1 
1 0 



1 1 
> = — [|A>±|B>] , i.e. | + > = — 

/2 /2 



and 



> " /2 



So 



H | + > = -A//2 


0 i 1 

1 0 




1 

1 


= - A//2 


1 

1 




l 











> 




1 \ 


0 1 




1 


1 0 




-1 






\ / 



> = -A//2 



(b) H in the {+,-} basis is -A 



A//2 



■ A | 



E = +A 



1 0 
0 -1 



Since H is diagonal, e 

T 

-BA 



-BH 



e BA 0 



0 e 



-BA 



and 



Q = Tre” BH = Tr 



BA „ 
e 0 



0 e 



= e BA + e BA as expected. 



A + > ; E = -A t 



Doing the trace with | A> and | B> gives 



Q = <A I 



0 ' 

0 

V J 



| A> + <B I 



e 6A 0 



0 e 



-BA 



B> 



|A> = [|+> + | ->]//2 
|B> = [|+> - | ->]//2 



1 r BA -BA, 

2 Ce + e ] 



1 r BA -BA, 

2 [e + e ] 



= i [ I 



(BA) 



n=0 



n! 



\ [e 6A ♦ e- BA ] 



I 

n=0 

1 0 
0 1 



( — BA) 
n! 



1 0 
0 1 



co n co 

]_ r V (BA) 

2 L L . n! n! 

n=0 n=0 



£ (~BA) 



0 1 
1 0 



1 r BA -BA, 
+2 [e - e ] 



r o A 
1 0 



cosh BA sinh BA 
sinh BA cosh BA 



These matrix elements are the same as those arrived at by rewriting 
in the {a,b} basis. Now performing the trace, 



e 6A 0 



0 e 



•BA 



Tr e BH = 2 coshfiA = e BA + e 6A 



In the end, the route to the partition function is unimportant, since 
the trace is invariant to a change of basis. 



(c) Using the { + , -} basis, 



(i) <m> = 1 Tr (me _BH ) = 1[<+ |m| +>e“ BA + <- |m | ->e +8A ] . 



Since 

<+|m|+> = ^ (<A|+<B| ) | m | (|A>+|B>) = 0 
and similarly <-|m|-> = 0 , therefore 



<m> = 0 . 



(ii) Similarly <+| |m| | +> 



-BA 



and <- | |m | | -> = e 



BA 



So < Jm | > = | u 



(iii) <(6m) 2 > = <m 2 > - <m> 2 = p 2 



H = - ia - pEa 
x z 



( d ) For E * 0 , 

(i) The energies then satisfy the equation 



-pE -A 
-A + pE| 



det (El - H) = 0 . 



So E 2 - (pE) 2 - A 2 = 0 



E = + 7(pE) 2 + A 2 



and Q = e 



p/(pE) 2 + A 2 + -b/(pE) 2 + A 2 



Then A(E) = -g 1 5,n[2 cosh(g7(pE) 2 + A 2 )] 



and the free energy of solvation is 
A(E) - A( 0 ) = -g' 1 ! £n[2 cosh(g/( pE) 2 + A 2 )] - dn[2 cosh(gA)] 



qtj 60 [ca + da ] 
(ii) e~ SH - l - z 
n=0 



n! 



where o = gpE and d = gA 



To evaluate the n— term, note 



2 2 2 

[co + do ] - c I + d I + cd ( a o + o o ) 

z x z x x z 



2 2 
(c + d )I 



So [eo z + do x ] n = (c 2 + d 2 ) n/2 I 



since o ,o = 0 
1 z x' 



for n even 



= (c 2 + d 2 7 n 1 ^ 2 [ca z + da x ] for n odd 



and 



BH _ v (7c 2 + d 2 ) n I 



e = I 
n=0 
even 



n! 



. v /(c 2 + d 2 ) n ( c d ■, 

l HI ° z + 

n=0 /2 2^ z fz 2 x 

n odd /c + d 7 c 2 + d 2 




| 

I 

n 



Bffssssgssssas 



3.22 



aligned with (or against) the field. 

As a final note, although <m> = - for small E, the fluctuation 

opti 

in m, <(6m) 2 >, is in general no longer equal to - unless m 

dBE 

commutes with H. Systematic analysis of the fluctuation in that case i s 
in fact equivalent to analyzing the quantum dynamics of the system, 
which is not the focus of this text. 



Bp - p/(1-bp) - Bap 2 , p = , V = L 3 . 



(a) For constant V = L 3 < ( 6p ) ^> ' / “Vp = <(6N) t > ' ' V<N> 



2 . 1 / 2 , 



, 2 . 1 / 2 , 



<(6N) > 



' 3<N> 



] = = V 

> a v u 



-3p 



3(By) J S,V k 3BvTB,V p'SBp'B.V 

= < N> (MP)" 1 

'■3p J B,V 



<N> 



(1-bp)' 



1 - 20ap( 1-bp) 



2 ’ 



<(6p) 2 > 1/2 1 



(1-bp) 



P ^ X-2Bap ( 1 -bp ) 2 

which vanishes as the volume becomes infinite. 



(b) At the critical point 

(ME) _ 1 

^ Jb " (i-b P ) 2 
Do some algebra! 



2gap = o and = 

3p 



2b 



( 1 -bp) ; 



2ga » 0 



p = -=7- and 0 
K c 3b c 



3b 



2ap o (1-bp c ) 2 2a (2/3) 2 



2Tb 

8a 



(c) Use the result of Part (a) with V = 100b, p = ^ and 



1 ilE 

x 8a 



where x will be close to 1 : 



, 2 . 1/2 

<(6p) > 



1 



Thus we have the table 



x 



< ( 6p ) 2 > 1 72 



1 .1 
1.001 
1 .00001 



0.38 

3-7 

37 



Lr 



(d) We're looking for "significant" density fluctuations in a volume 
V = - (1000A) 3 . Typical molecular volumes are - ( 5 A) 3 = b 
¥ = (200) 3 b = 8 x 10 6 b . 

< ( , 21 /2 

To estimate significant - — ^ — , we checked the CRC for the 

index of refraction of distinguishable (by sight) fluids: 



H 2 0 = 1.33 

n-Hexane =1.37 or ~ = 3% . 

For small differences in p, we can consider nap. So take 
<(6 ) 2 > 1/2 

significant = 0.03 • Solving for the difference in 

B = B Q /x at p = p c = , we get: 



- 10-3 -I . 

/6 P 

c , tx , 21/2 -1 

Then x = 1 + 10 5 f SlRI 1 = 1 + 10 “3 

P 



In other words, critical fluctuations will be easily observed 
optically if you are within 0.1? of the critical temperature. 



A * B , N A + N B = N . 

We can make an analogy with the ideal gas analysis of Sec. 3.6. 

<{ n a - <n a » 2 > = <n 2 > - <n a > 2 

N 

= I ^- <n ai n aj > " <n ai ><n aj > -' ’ n ai = 1 ’ if raolecule i is in state A 



= 0 , otherwise 






zero since different molecules 
are uncorrelated 




= N<n >(1 - <n a1 >) , since each molecule is equivalent 

“ N<n ai X1 " n a 1 > = Vb N ' 
where 

X A = <N a >/N , and x B = 1 - x ft . 



Chapter 4. Non-Interacting (Ideal) Systems 



5 Begin with the equation at the top of page 9*1, 

BA = -inQ = / du g ( u) ) in [e SMu/2 - e ' 6K “ /2 ] 

0 



2 w 0 

- (^-) / dw Co 0 ’ 1 in [e SM “ /2 - e' 6 ^ 72 ] . 

“0 0 

As long as 6^/2 >> 1 , we can approximate this integral in the limit 
of B + “ as 

2 “o 

lim BA = (— E— ] j du u D 1 BMu/2 
T+0 ui„ 0 

(8 — ) ° 

ND u Q M 

= 2( D + 1 ) 8 ' 



_ 3inQ , ND “o M 

3 ( -6 ) ' N , V 2 ( D+ 1 ) 

Similarly , 



C = (11) = - k B 2 i_ 1 

V l 3T J V V 36 3 (-6 ) \ 



2 3 r" a 6Hm/2 x -ghai/2 

k B 8 3g C d “ - BM|fl/2 - e - gMu)/2 



Mu/2 ] 



2 “o 

■ - v 2 f ' < 0 * ■” <“-> 2 . W 

Change variables to x = gKu and let k D 0 n s Mu : 

B D 0 

c = k Jill M f° D/ ' £1 H x 2 e X 
V %° 9D 0 8 « g DX IV 

0 /T 

p T n , D D+1 X 

■ V ° ( fc ) / hrrz “ • 



In the limit T + 0 , we encounter the integral 



“ D+1 x <o 

f x e D+1 -1 

J — p dx = * 

0 (e -I ) 2 e -1 0 



(D+1) / 



0 e X -1 



= 0 + (D+1 ) r(D+1) 



l( D+ 1 ) 



(D+1)! C (D+1 ) ■ 



Gamma Reimann 
Function Zeta Function 



For D = 3. we note C( 1 *) = 11 (9° > and 



12 „„ fL-1 3 

E- ir Nk 0 (, Q J 



D 



<n.n.> = 
i J 



1 a 2 s 

H 3Be.38£j 



3(-&£j) 



<n. > <n .> i 56 j . 
i J 



= <n.> 

l 



i = J 



(5 <n.>) 



since n. and n. are uncorrelated 



since <n^> = <n^> 



Therefore, 



j = <n. ><n .> ( 1-6 . • ) • 
4 ij i J U 



, n The expression for the electronic partition function of an H atom, 

B£ 0 B£ 0 /If e BE ° /n + ••• is oniy exact in the 

q int (T) “ g i e + g 2 e + -'- +g n e 

case of an isolated atom. When the atom is put in a box or if there are 
other atoms in the universe, the convergence of energy levels 
->E = 0 asn + “>is lifted. The reason for this is that as n increases, 
so does the energy and spatial extent of the atom. The expectation 
value for the repulsive energy of two atoms then has no upper bound, 
will increase as the spatial overlap increases. 



.15 For an ideal mixture 



V 



V 



Further 



38-nQ -, 

3N . p , V , N p 



-BP, 



V 31 



Using 8 ,nN! = NlnN - N we get 



14 . 1 6 



-S.nQ = - N A 8 ,n q fl + !ln N fl ! - N g In q 0 + in N g ! , 
so that 



P A = u B implies -In q A + In H ^ 



-Zn q g + In N 0 



or 

N. 



(a) 



A -BAe 
— e 



B 4 B 



„ 1 N 1 , , N 

Q N! Q N! q A + q B 



1 v N ! N A (N “V 

N! N =0 N A !(N_N A )! ^ ^ 



n a n b 



L N !N ! q A q B 

n a ,n b VV 



l exp[-BA(N A , N b )] 



such that 



N A + N B = N 



/ l \ , . 3A 3 A 

(b) p A - p B implies 3 — = 3 — . 



<N a > + <N 0 > = N = const implies ^ . 



0 ) 

( 2 ) 



Equations (1) and (2) imply 

3A 3A 3A 

3<N a > ' 3<N a > ” 3<Ng> • 

Therefore, quite generally, A(N A ,N B ), is minimized when 
N A + Ng = N are partitioned to satisfy chemical equilibrium, 
specifically, with the ideal expression in Part (a), we have 

0 = gf^ I" J.n[N A !(N-N A )!] + N a Jin q A + (N-N A Hn q B ) 

= - S,n[N A /(N-N A >] + S,n(q A /q B ) , applying Stirling's 

approximation 



4.17 



or associating the solution to this equation with 



4.18 For an ideal gas of structureless particles 



<N„ > = N - <N > , we obtain 
A B 



<N A >/<N B > = V q B • 



Q = I 



V N B 

(n a + n b = n) 



n a n b 

Q ft Qd m 

_JL. = (n + q ) VN! as shown above . 

N ft ! N d ! V4 A 4 B' 

A D 



n a n b 



'-ft 



<N A> ■ Q 



N A’ N B 
(N a ♦ N b - N) 



V 1 N B 



A N a ! N b ! 



n a' n b 

(n a + \ - N) 



A N A ! 






2 2 

W k 

"k 2m 



’ - = L ( V + + n A 



(a) Then assuming a macroscopic volume, 



8(u-e k ) 

BpV = Ins = l Jtn[ 1 + e ~ ] = 
k 



2 , 2 , 



(air) 0 0 



/ WWl + e 5( ^ M k /2m ) 



]dk 



4 V , m -, 3/2 2 , -x 2 

2 l " " o J J x £n(1 + ze ) dx 
tt 2Bh 0 



So 



- 73 / 



A r 0 

1 



? - c 
X tod + ze x ) dx 



3 f 5 /2 (z) • 



z = e 



Bu 



, = r 2ir B M 1 1 /2 



C 3 !,nQ'| 

= q A 1 3 q A J Q B .N 



Substituting the series for S,n( 1 +x) , 



N 



, , N « 

, a (q. + q R ) , 

q a Un ( Hi )] q B - - ■ m b - 



J q n ,N " q A (q. + qj 



Since A and B play symmetrical roles, <N g > = q g — +~q“) ‘ 

A B 

Hence, 



<n a >/< n b > - q A /q B 



Similarly, 



n a. n b 

2 q A q B 



1 3Q ^2 



( — ) 



«»«- <V> 2> ■ <k » 2> - <n a >2 ■ 5 , l . "* w 7 [ 

A B 

(N a +N b = n) 



q A f 3 Q a 2 n A A_ 1 3 Q_ _ ^A f 3 Q 

Q 3 q. '"A 3 q/q n ,N ' n 2 1 3 q„ j q D , N ' q A Q o 



*a a ( 9Q \ 

~ A 1 A A^b 



4 a m b’ 



3q ft 2 Q 3q A Q 2 A 



. 0 !^* 1 | = 0 [ N - q » 

q A 3 q ft Q 3 q ft q g >N A 3 q ft q B> N A + Q B A (q fl + q B > C 



q A N C V Cq A + V ] = <N a ><N B >/N 



K 7 



= 7 1 ^> K+1 f I 

V it k =1 1 

- I (-1) k+1 z k / k 5/2 . 
k=1 



To calculate the density, 



<N> = I <n, > = Y — 7 v 

k k k ^ 



2 -kx , 
x e dx 



+ 1 



(2 tt) j 0 



/ 4irk 2 



2 2 

-MkV2m 



1 + ze 



2 2 

-H k /2m 



dk 



V 4 



■i / 



2 ze 



A /tt 0 -x c 

1 + ze 



dx . 



So 



pA 3 = Ji y (-1 ) k+1 z k f x 2 e' kx " dx 
Sir k=1 0 



- I (-1) k+1 z k / k 3/2 
k=1 



f 3/2 (z) . 





So pA^ = (S,nz )^ / ' 2 



H/3/H 



! 



-B(e-e F > 



BE- 2 2 2/3 

and z = e * where e = (M /2m)(6ir p) 



Be p -x 



(b) BpA 3 = r U) = — j dx x 2 in(1 + e ) 

5/2 A 0 



= e F J Q d e' f3/2( z > 

To get the correction term in (BCp) > we calculate 

p 

order t : 



= 2p e 
s 



Similarly , 
2 



n 



( 2ir ) ~ 
2 



f dk [1 



B(e,+ £„) 



k F\-1 2 



( 2-n ) “ 



f dk 



” 6(E k +£ F ) 



"BA 



4 it J dk k 2 exp(-f5K 2 k 2 /2m e ) 



I dt 4 ^ 2 

(e Z + 1) 



4.2 t 

t e 



{ dt t 2 

(e 6 + 1 r 



0 since the integrand is odd, and 



So SPA 3 = e„ j dp',— [(B'e p ) 3/2 + | (B' EF ) 

0 3A 



- 1/2 



+ ..•] 



4 (Be F ) 3/2 [ | (Be F ) + \ <Be F )“ + ••• 3 



3A 



Hence p = f- p £^[1 + (B£p) + •••] 



4.21 (a) Simply plug in the Fermi di 



-gU-Ep) 

(b) Since e = x << 1 , 



then — = 1 - x 

1 +x 



distribution given on page 96. 



So, 



p = 2p s - 2 p s C 1 + e 



-B(e-e ) _ 1 r . ^^ e-£ F ) •>-, 

F ] = 2p [1 - (1 - e J] 



j 

( 

* 



(2ir) 

-Be,, . m 



- e F 2 P4) 3/2 

2itBM 



" Be F 3 
= e 2/A 3 . 



Therefore , 
pn = p g /A 3 4e Be 



5 £ “f 1 - 

The pressure is not zero at T=0 because the exclusion principle 
requires all but one fermion to have non-zero momentum. 



4.22 <N> ad = j 



<n.> = p V 
J s 



-B(-e-u) = e By Ve Bc 

K s 



adsorbed 



degeneracy classical 
limit 



ad 



^ld . gBu e BE 
V H s 



Also By = Hnp g A 3 for the classical gas, and y has the same value as in 
the above equation for p gd . 




Hence , 



p = e / a 



which implies 



/ , 3„Be 

p ,/p = p A e 

p ad s 



4.23 (a) E = E A + + E c 



then 



Q = X e 

ABC 



•B(E a+ Eg + E c ) 



-BE. -BE -BE 

Ale A ) (X e B ) (X e 

ABC 



> = 9„ 9 n 9, 



A 4 B 4 C ’ 



e 



e 



3S,nQ Wnq A . 3toq B , 3toq C 

3(-B) 3(-B) + 3(-B) 3 <-B) 



, and 



3 2 S!,nq 3 2 S.nq 3 2 8,nq c 

C , = -k Ba<E>/3B = k B^[ — p— + p - + T ] 

v B B 3B 2 3B 3B 



r (A) . -(B) . -(C) 

C V + c v c v 






E n = zero point energy 



-B(E a+ E 0 ) -6E 0 -BE a 

3 = e i e 



PC*. P*-*A 

35,n Q 3 [ _ BE + ln l e A ] = E„ + l E e /Q , and 



0 A A 



C V = ' k B 6 *■ 3B 



2 &£) - - M 2 1 * d ^ • 



B p 3B ; L h “A 



which is independent of Eg. 



In other words, fluctuations produce C v , and Eg has nothing to do with 
fluctuations . 



-Be, -Be 2 

(b) q elec ‘ g 0 + g 1 S + g 2 e 

where energies e^ and e^ a he relative to Eg , our choice for the 
zero of energy. 



(eiec ) k x 2 3_ 

C V B ag 2 160 q elec 



= k T 2 { — 
B 1 q 



2 2 
[g^e + S 2 e 2 e ] 



pE 1 B 2,2 i 

[g^e + g 2 e 2 e ] t 



Trans 4 elec 4 nuc 4 rot 4 vib AB 



(trans ) (elec) _(nuc) _(rot) p (vib) 
C V + C V C V V V 



Now 



cJ ,leo1 . - &> 2 ] 



2 -Be, -Be 

2 2 , g i e i e + ••• S^e.e + 

<e > - <e > 2 = -LJ _ ) _ fJJ 

g 0 + 1 g 0 + "- 



7 ,-50,000/300 -100 

' ® N 0 



(elec) 2 -100 

^ k B B 6 at roonl temperature, which is negligible 



Because of the high energies of excited electronic states, the 
contribution to C y of excited states is negligible. Neglecting 
the excited states completely is equivalent to assuming 

( ©1 6C ) 

C V = 0 , which is therefore a good approximation. 

Also, the above calculation shows that the ground state 
degeneracy doesn't enter into the heat capacity. 



<E> - A 



q elec = g 0 e + g 1 e 



. -100 " pE 0 
< e e 



Ignoring these terms doesn't affect q elec 
significantly and also won't affect either 
<E> « 3&nq/3(-8) or A = -ilnq/B . 

On the other hand, g Q will affect S because even 
though the average energy doesn't care about the 

degeneracy, <E> = g 0 E 0 e B£ 0/g 0 e~ Be 0 „ e _ fche free 

energy is affected. 



" B g 0 + E 0 • 

V + k B T lng 0 



= k B tag 0 ■ 



Thus the degeneracy accounts for an additive term to the entropy 



1) 2 4 Q = q N q N /N ! 

w 4 trans ■ int 



ft 



If we write 



<5 int = s o e (2I A + 1)(2I B + ^ q rot q vib 



‘Vot ’ T/0 rot 

BHiOq/2 -6Hu 0 /2 _ 1 

q .. = [e - e ] . 

^vib 



we are assuming: 



( 1 ) The volume is large enough that q int is independent of volume. 
This is important when distinguishing C p from C v , since the 
difference will be that from q trans . i- e -- the ideal gas value ' 



(2) Excited electronic (or nuclear ) states are sufficiently high in 
energy to be ignored in Q. 



( 3 ) Internal nuclear motion is well approximated by a rigid rotor- 
harmonic oscillator, which neglects: 

- anharmonicity in the vibrations 

- centrifugal distortion 

- rotation-vibrational coupling. 



(4) The high temperature formula for q rot assumes that 

0 << T (i.e., 12.1 K << 300 K) . 

rot 



From the partition functions, we obtain 



<E> = -3S,nQ/3g - Ne Q + NkgT + N -jr coth( BKu> 0 /2) + 3/2 NkgT 



Therefore, 



0 p p o 

C v = 5/2 Nkg + n(~) cscfi ( BMiOq/ 2)/ kgT , 



and C = C, . + Nk„ 
p V B 



= 7/2 Nk B + Nk B (2ipjr) 2 C3Ch 2 (B»ui 0 /2) . 
B 



Furthermore, 



<E> - A 



Nk { 5/2 + — eoth(Bhoi) /2) + {.ng + S,n(T/0 . ) 

d d u nuc rot 



- Jtn[e 



gKm n /2 -gMm n /2 



2ir(m + m ) 

] + tn[ ^ ( V-- 3 3 

h 3 g 



S,nN + 1 } . 



The necessary numerical input is: 



T = 298 K 6 hu, 0 - 



g nuc = (2 r'> (2 I + D 



1 = JL = 4.06 x 10 2 cm 3 

N p 



Br 81 , Br 79 



yielding 



S = 51.7 cal/mol K , (Note: CRC value of 47.5 neglects the in g 

contribution of 4.14 cal/mol K) n 



Cp = 6.97 cal/mol K 



4.26 Q( N , V, T) - 3_ 



V_ 1 2irmi 3/2 
n v- a ) 



V q N „BpN “ (qe 6w ) N (qe gu ) zV 
^ mi e L mi — e e 






e Bu f 2irmi3/2 

Z ' 7 l T> 



<N> = ^ n — = V z implies z = p 

3Bu 



Hence, p = p/3 , the ideal gas law. 



/<(^p) > = 



/<(Sn) > 

<N> 



<(^) 2 > - Hr = v fe = Vz = <N> 



/<(jp) > _ i 



k n T/pV = 2.01 X 10 



(d) To second order in 6N 



>(N) = £nP(<N>) + 1 («JN ) 2 ( 3 iln P ( - ” i-) 



3N N = <N> 



i.e., assuming a Gaussian distribution around <N> with variance 
<(^N)> 2 , we get 



P ( <N> ) 



e 10 <N>/2 _ e 10 j which is very unlikely 



Chapter 5. Statistical Mechanical Theory of Phase Transitions 



5.1 Since the number of nearest neighbors in a D-dimensional cubic lattice 



is 2D, 



E 0 = " JN ( ^ • 2D) = - DNJ , 



where the V 2 corrects for double counting pairs. Note that for a two- 



dimensional triangular lattice, each spin has 6 nearest neighbors, so 



that in that case 



S 0 = ~ JN ( \ 6) = " 3NJ . 



N 

- J l 



s . s . „ 

1 1+1 



N 

Jib. 

i = 1 



where b. - s^.^ = ± 1 is the bond between spin i and i + 1 . The b.'s 



completely determine the Sj’s up to the value of s, . This degeneracy is 



related to the fact that only N-1 of the b^'s are independent: 



n b. = n sf - i 
i=1 i=l 1 



However, as N ■* ® , the contribution in energy due to the nUI bond 



N-1 n-1 

b = n b is small compared to the other bonds 7 b So 

i=1 1 i 



Q - , I 



BJ I 3.S, „ 
i=1 1 1+1 



BJ I b. 



ls.= ± 1 ] 



“ , I e 

N-+® |b_. = ± 1 } 



- n ( I e gJb ) - (e BJ + e‘ SJ ) N - [2cosh( B J)] N . 
j=1 b — 1 , 1 



Note: With some patience, the reader can compute Q for finite N without 



neglecting edge effects. Also of interest is the transfer matrix method 



illustrated in the solution to Exercise 5.21. 



I (2n -1 ) - J l' (2n i -1)(2n j 
L-1 i.J 

N _t 

„H l n. + tnNH -HJ I n.n + - 
i-1 i.O 



N 

same as z 1 



esult , the answer is 



(B.N,H;J,m) = exp[BN(mH - Jz 



U = 2mH - 2Jz, e = 4J and z is 
to any given lattice site. 



g 2(ByH + B 

tanh (ByH + BzJm) = ?( r u H + 1 







In this figure, only positive m are plotted since the curve is symmetric 

(MF) (MF) 

about the temperature axis. Note that kgT/2DJ is T/T c , where T c 
is the mean field theory prediction of T c . 

5.8 In the mean field theory, different spins are uncorrelated. Thus 
(£)mf ~ ^ ~ 2 *^0' 

i ij 

= X^( s *)mf ” 2 X/ > 

i i.J 

The desired result follows from the association {si) UF = m. Al- 
ternatively, you may begin from the perspective of Sec.(5.5) where 
the partition function in this approximation is shown to be 

InQ fz IhQmf — ^(A-E^p 

where —/3(AE) UF is given by Eq.(d) of p. 138, and 

Qmf = [2 cosh (fifiH + /3zJm)] N . 

The desired result for the internal energy is obtained from differ- 
entiating this approximation for InQ with respect to —/?. Here, 
in the differentiation, note that m is a function of /3; in fact, from 
Eq.(b) on p. 136 with AH = zJm/p, we find 

~~ = ( nH + zJM)( 1 - m 2 )/[l - PzJ{ 1 - m 2 )] 

for T < T c . For T > T c and H = 0, mean field theory gives 
m = 0, and thus predicts zero internal energy in that case. This 
prediction is incorrect since (although the average coupling to the 
field is zero — i.e., there is no broken symmetry' and long ranged 
order) there are still correlations between nearest neighbor spins. 



18 g(K ' ) = 2g(K) - h(K) 



where h(K) is analytic at K = K 



K = K (K) = 3/8 £n cosh(4K) 



K CK c ) = K o , 
and h(K c ) = g(K 0 ) 



Assume C « |K - K I “ near K = K , 
1 c 1 c 



that is, C = ( a less singular part } + a|K - K Q | _a . 



The second term with a > 0 dominates as K -> K 

c 



Then since c = g(K) , 

dk 2 



we can write g(K) = { a less singular part } + b | K - K c | 2_a 

where b is a constant. Plugging this into the first equations above, 



(a less singular part} + b | K ’ - K 0 | 2 “ = 2[{a less singular part} 



+ b|K - Kj 2 “] - h(K) . 



As K -> K q , the dominant contributions are 



b|K K„ | 2_ “ = 2b | K - K„ 



Taylor expanding k’(K) about K = K q gives 



K ■ K « * Sk ■ « - V 

c 



| § | Kq . (K - K c ) | 2 -“ - 2 1 K - K c | : 



9K | | 2-a 

8K 'K 1 



2 , which implies (2-a) J,n 4r-L = In 2 , 



or 



5.21 



sic 

a = 2 - An 2/ fcn |^-| K 

c 

Note that the actual value of h(k) is irrelevant. 



(a) Tr(q N ) = l 



(b) Tr (q N ) = Tr[ 





N 


(S j + 

e J 




n 


*2 * ” ’ 


S N =±1 J=1 


N 




exp[ l (hs 


3^, . . , 


■ ,s =±1 




2 


N 


x o] 


N 


t N X 

1 a" 0 


+ 

0 X 


] = Tr [ 


In 

1° 



J J + 1 



)] = Q 



] = X” ♦ x" 



(o) q 



a c 
o b 



| XI - q| = ( X-a) ( X-b) - c 2 = X 2 - (a+b)X + ab - c 2 = 0 



implies 



X = ± /(^) 2 ♦ c 2 



Thus , 



\ = e^eosh(h) ± J e 2 ^sinh 2 (h) + e 



e K [cosh(h) ± J sinh 2 (h) + e ] 



-UK 



Further , 



X N 



Q = \ + + X_ = x"[1 + (y 1 ) ] , so 



8,nQ 



X N 



JlnX 



+ + i ta[1 + (— ) 3 ■* * n *+ > 



as N 



Hence , 



j 2. 

= Jtn ( X ) = K + S>n[cosh(h) + /sinh (h) + e ] 




sinh(h) 



sinh (h) + e 



Thus 



Q(K,h,N) = 



(h/2)[s +s ]+h+K[s +s ] (h/2) [s +s ]-[h+K(s +s ) ] 

l (e 13 1 3 ♦ e 13 1 3 ] x 

(odd s^ ) 

(h/2)[s N _ 2+ s N >h+K[s N _ 2+ s N ] (h/2)[s N _ 2 + s N ]-[h + K(s N _ 2+ s N 

. . . x [e + e 



= [f (K,h) ] N/2 Q(K' ,h' , N/2 ) . 



(h/2)(s +s )+h+K(s + s ) (h/2)(s 1 +s ) - [h+K(s^ +s )] 

e ^ J + e 3 ' '' 



(h'/2)(s +s ) + K's s 
= f(K,h)[e 1 3 13] 



f(K,h) = 2[cosh 2 (h) cosh(2K+h )cosh(-2K+h ) j 1 ^ . 



sinh(h)cosh(h) 



cosh(h) + /sinh 2 (h) + e 



lim <s > = lim sinh(h) 
h->0 1 h+0 



1 + cosh(h)//sinh 2 (h) + e ^ 

/ 2 -4K 

cosh(h) + /sinh (h) + e 



By multiplying Eqs (1),(3) with Eq(2) twice quickly yields 



the possible states of s^ and s^: 



, h'+ K' 2(h+K) -2K 

s 1 - s 3 - 1 : f(K,h) e = e + e 







(c) g(h,K) = 1 S,nQ 

g(h,K=0) = i 8,n(e h + e h ) N 
N 

= in(2cosh(h) ) 

As written in part (a) and (b), flow from g(h,K = 0.01) will be 
towards the h-axis, away from (h,K) = (1,1). V/e would like then to 
invert the parameter flow, i.e., find h = h(h',K') ; 

K = K(h ' , k' ) . 

Using these inverses, our errors will diminish with each iteration, 
since then g(h,K) = 1/2 g(h',K') + 1/2 f(h',K') . 

One method to get reversed flows (from low K to high K) is to first 
flow forward f, om (1,1) (to lower K values), and then retrace the 
steps taken back up to (1,1): 




h 



On the reverse trip, d-e-f, we can use our low K "accurate" 
approximation for g, the error halving with each iteration. 






g 


h 


K 


2.00285 


1 .00 


1 .00 


2. 41 


1.94 


0.469 


2.84 


2.81 


2.23 x 10 -2 


2.86 


2.86 


7.09 x 10 -6 


2.86 


2.86 


7.5 x 10~ 12 



Generate starting at (h,K) = (1,1) 
and moving down according to flows. 



Start at bottom with 

g(h = 2.86, K = 7.5 x 10’ 12 ) = g(2.86, 0) 

=u Sln[2cosh(2. 86) ] , 

and iterate upwards using 



;(h,K) = 1 g(h' , K') + j Hn f(h,K) . 



The convergence to h = 2.86, K = 0 exhibits a stable fixed point, 
and in fact any h = h Q , K = 0 is a fixed point in flow. 



The exact answer is (from Exercise 5.22) 



i>nQ(K-1, h-1 _ ) _ = K + j,n(cosh(h) + /sinh 2 (h) + e"^ ) = 2.00285 



which matches the RG value of 2.00285 above. 



An analytic inverse in closed form can be found by using e^ 1 and 
e" k as variables, which reduces the problem to the finding of the 
roots of a 4^5. degree polynomial. 




24 Parts (a) and (b): The rules for construction of legal configurations 



admit only three energy levels 



if any partition has other than two atoms. Now, given that 



there are two atoms per partition, we can have 



Having two atoms on the same side of a partition means that the 



two atoms in the next cell must belong to the next partition 



The degeneracy is 2, and the energy is E = N-O = 0 



(iii) The only other noninfinite energy states are those where each 



partition has an atom on each side 



mea ns that the next partition is either 



and the choice between these two possibilities for each partition 



is unrelated to other partitions, e.g 



These states have energy E = Ne, and the degeneracy is 2 ] 



As a result of the preceding analysis 



Therefore 






Hence 






<E> 


1 3£nQ 


Q “Be 
e , 2e 


<N> 


<n> a (-b ) 


1 0, 2e _Be 



-"<V/ |Sg>, 



(e) A phase transformation is associated with g - — in 2, or 



U J dE e BE /2a 2sinh(g»4 2 + p 2 E 2 )-pE//a 2 + p 2 E 2 



‘0 " k Jln2 • 

D 



This problem illustrates how phase transitions can arise due to 
competition between entropic and energetic effects. 



/ dE e _BE /2 ° 2cosh(/i 2 + y 2 E 2 ) 



As 8 + 00 , coshBx, sinh8x -> e^ x 



e , so 



U / dE exp[-0E 2 /2o + 8 /h 2 + p 2 E 2 ](pE //t\ 2 + p 2 E 2 ) 



5.26 (a) Q = (eA) P/2 l e*P[ l * u u ]/ dE exp[-BE /2o + (ep I u.)E] 

P "“ lu.l i-1 1 — i=1 



/ dE exp[-6E 2 /2a + b/a 2 + p 2 E 2 ] 



The integral can be evaluated to be 



>{0/2a[f (ep l u )] 2 ) / dE 



-B/20[E-o/B(eu I u.f 



P / dE e BW E (yE//A 2 + p 2 E 2 ) 



/ dE 



- exp[ ( \ u )( I u )] • 

2P 2 i = 1 j=1 J 



l < u lUi+1 + (#M Vj 
_ /2™ lim {(eA) P/2 i e 1= ' 2P 1>J = 1 ] 

g P->» { u . } 

l x J 



where the last equality defines W(E) and the weighting function 

-bw(e) 

& . As 8'*“ , the weighting function becomes a delta function 

K 6(E - E m ^ n ) where E^ n is the E that minimizes W(E). Exploiting 
this idea is a method of "steepest descent." Note 



3W = E 
3E ° o 



/2 2l2 

/A + p E 



which implies E . = 0 or op 2 = /a 2 + p 2 E 2 . 



(b) From the text, we know 



l [< u.u i+1 + hu.] 

lim V / ai P/, 2 i = 1 
Q(E) = .1 e 



In the latter case, 



22 . 2,2 



n = o u - A /p which has the real solutions 



= 2 cosh[B/A 2 + p 2 E 2 ] 



E ± = ± 



A_i 1 /2 
2 
V 



Hence , 



when op /A > 1 . When this condition is met, we have broken 

symmetry. In other words, E + starts contributing to W(E) at 
2 2 

op /A - 1 . For 0 p /A < 1 , <m> = 0 , while for op 2 /A > 1 , 



g->“> , <m> * 0 . 



Thus, we're looking for bistability in A like this: 




Thus 



This emergence of broken symmetry will cause 
<(6m) 2 > to diverge at o = a ^ = ’ More is saifi about this 

broken symmetry in the isomorphic problem, Exercise 5.27. A 



5.27 (a) Doing the sum over spins with the transfer matrix method (Exercise 

5.21 ): 



J exp[ T Bhs. + gJs.s. .] = Tr(q N ) , where q - ( _ 

, ,1 ji 1 ii+i - „ 



Nov;, the eigenvalues of M = 



e Bh+gJ e ~BJ 
e -BJ e _ Bh+BJ 



where A > A_ are the eigenvalues of g , and 



Since A is even and A ■* “ , A will have a bistability in h when 
h->» 

9A I 

9hlh' ^ ® f° r ^ 0 near zero; i.e., A is concave downward at h = 

0 when 3A/3h' < 0 for h' near zero. Hence, 3A/3h = 0 can be used 
to determine the critical temperature. Alternatively, note that 
taking the B- >c ° limit gives 



! (h; B > N) = e 



-BNh /2o r ,N .N 



lim cosh(gh) = e e ! h l , lim sinh(Bh) = sign(h) -e B l h I , and hence 
S - *-” B^“> 



For large N, therefore, 



= N[^--J--jS.n (cosh(Bh) + «4inh 2 (6h) + e 6 )] 
2o 3 



(b) At this point, notice that A(h;B,N) is even in h , i.e., if h 
minimizes A , so does -h . Further, 



and does show bistability. However, we don't get an equation for 
B crit from this analysis. 



= N [ — - sinh(Bh) e 2 ® J ] 
- + n 



lim 



since sinh(Bh) = Bh . 



This derivative changes in sign as a function of 3 when 



B . . e 
H crit 



2B crit J 1 



which is an invertible function in 6 since it's montonically 
increasing. 



Another way to view the analysis of this Exercise is to solve for 

M - o to find the stable minima in A(h) . This equation can be 
8h 

written A'(h min ,8) = 0 , an implicit function for 3 = e( h min ) • 

g is then the limiting value of 8 = B(h j* ) . as h mln + 0 > 
crit 1 



and so we look for - 0 and solve for B crit - 

d 'h->0 



Chapter 6. Monte Carlo Methods in Statistical Mechanics 




To perform Exercises 6.10, 6.11 and 6.14-6.16, you will collect 
tatistics from long Monte Carlo runs. These exercises help illustrate what 
must be done to obtain reliable information from such calculations. In this 
egard, efficient and fast computer codes are important, and as emphasized in 
he text, much can be done to improve the programs listed in the text. At the 
ery least, one should use a compiled version of the codes. As a standard of 
xcellence, however, you might want to view Jeffrey Fox's Monte Carlo of the 
Ising magnet.* It utilizes binary arithmetic only and machine language, and 
it runs from 10 2 to 10 ^ times faster on a PC than the Ising model program 
listed in the text. As well as providing rapid access to good statistics, the 
^reater speed permits one to view rather large systems in real time. The 
iews are both instructive and pleasing to the eye. 

Some new books on Monte Carlo methods in particular and computer simula- 
ions in general are worthy additions to the bibliography of this chapter: 



M. H. Kalos and P. A. Whitlock, Monte Carlo Methods Vol. I: Basies , 

(Wiley Interscience, New York, 1986); 



H. Gould and J. Tobochnik, An Introduction to Computer Simulation 
Methods, Parts 1 and 2, (Addison-Wesley , Reading, 1988); 



M. P. Allen and D. J. Tildesley, Computer Simu l ation of Liquids , (Oxford 
U. Press, Oxford, 1987). 



The last of these is especially pertinent to the material covered in Chapters 
7 and 8. 



Fox, "Fox's Ising Model Simulation for PC Compatibles," placed in the 
ublic domain on 3/25/88. 





6.10 Above T 



decays more or less exponentially as i gets 



farther from j , and random thermal fluctuations destroy long ranges 



correlations. In fact c. , - 0 for i far from j until T - T„ from above 



at which point * 0 



Below T c , we have broken symmetry, and as T -> 0 , we find 



is significantly greater than zero. This is an indication that 



we are frozen into either the up state or the down state. One result 



this broken symmetry is that as our runs apparently show <s,> •* 1 



and hence <s : s,> + 1 as well, we will observe c.. ■* 0 as T + 0 



Spontaneous 



magnetization 



Correlations in the 



Ising Magnet 



Representative results for exercise 6.10. The lattice was equili 
for 100 passes, and then statistics were compiled for 5,000 passes 
Spin-spin correlations were measured up to the fifth nearest neig 




6.11 The effect of the interface is to introduce long-ranged correlations 
between spins along the interface, even while T << T c and the 
correlations in the bulk are only short-ranged. However, as defined in 
the code given in the book, the ground state of the 20 x 20 lattice, 
even with the biasing field on, has all the spins aligned the same way 
(with the exception of the biased spins). Simply setting the initial 
configuration to be an interface would sample a metastable well. To 
observe the interface fluctuations we can increase the bulk of the 
lattice relative to the interface by changing the dimensions to be 
40 * 20 . 



Note that in implementing the calculation, there are four columns 

of spins each that oan be considered bulk or interface. In our 

calculations the dependence of <s> and <SjSj> on the distance from the 

biasing field Is averaged over. Correlations shown in the graph are 

<s s.> - <s . ><s . > . 

1 J i J 




Neighbor distance 



Fig. 



These representative results were obtained on a 40 < 20 lattice at 
k B T/J = 1-2 usin g 100 passes to equilibrate and 10,000 passes to 
measure correlations. The biasing field keeps spins 1-20 up and 21 -A0 
down in the top row. Spins in columns 1, 20, 21 and A0 are considered 
to be in the interface, and columns 10, 11, 30 and 31 were taken for 
bulk measurements. The upper line shows correlations along the 
interface, and the lower line shows the correlation of a spin at the 
interface into the bulk. Correlations between bulk spins are - 10 - ^ o: 
10 and are not discernible on this graph. 





6.12 A straightforward algorithm is to choose a large enough, or 

equivalently, to scale the range of the uniform random numbers [0,1] t 

[ - x m „,,, x m = v ] so that most of the Gaussian will fit in this range, 
max in a x 

trial x's will be accepted with a probability <* p(x) , normalized so 
that for all x, prob(x) < 1 . The trouble is that the fraction of 
accepted trial x's tends to be small. This fraction is just that of the 
area under the curve to the rectangle: 




In other words, if we pick x uniformly from the interval, we will the 
reject most choices. This inefficient method does not exploit any for 
of importance sampling. There are a couple of possibilities for 
improvement . 

2 

One is to imagine that the distribution p(x) = /a/n e is due 

2 

to an energy E(x), where gE(x) = ax , and to employ the usual 

Metropolis algorithm with diffusive motion. Then trial samples will 

tend to come from the regions of lower E(x) and higher probability. The 

drawback is that one random number x^ + -| is correlated with the previous 

one x^ from which it made its move to x^ = x^ + Sx , with 6x 

small . Nonetheless , the calculated moments should be correct . 

Another method, equivalent to the "force bias" Monte Carlo of 

problem 6 . 13 , is to change it , to reflect the more likely regions 

vv 

p(x), thereby increasing the A , 's . 




The first algorithm described above corresponds to w 



constant 



As we said, that scheme is not very efficient, if 



instead we make a change of variables, such as x = f(y), with y picked 



uniformly from [0,1] then the probability of picking x in [x n ,x. 



where y - f '(x). So if we want to achieve 100? acceptance 



1 , then our trials must be chosen with 



and thus 



which implies f 



f is the inverse of the error function. Note this scheme only 



generates positive X; the sign should be assigned randomly 



In some problems, such a perfect solution is not possible. But as 



long as A yy , increases and begins to approach unity, we have made good 



progress. In the present context, any change of variables favoring 



sampling smaller x is better. For example, let 



where y e [0,1] is uniformly picked. This implies 





So the x are picked with probability 



ir « , and so we must accept with probability 

x . ,2 



A cc 
X C 



(1 + cx) 



(1 + cx) 



/FTa 



equivalent to any other pair. With detailed balance, designating 
E i > E . v' = v+1 , we therefore have, 



v v 



-6AE , 
= e vv 



-gMoo 



1-p 



to preserve w x . The proportionalities are chosen so that n < 1 and 
A v < 1 . Again only positive x are generated and the signs should be 
randomly chosen. Note only one index is retained because any given 
choice is completely independent of previous choices. 






P = 



1 + x" 1 1 + 



Note: in the limits Hu + 0 , p 1 /2 , 

and hoj -» 0 ° , p -* o 



6.13 (a) The condition of detailed balance, w ,/ w , = exp(~3AE ,) , 

vv v v vv 

with AE t =E,-E , gives 
vv v v 



6. 1H (a) Q = / dE[Tr 



?(h - m(E + E ) ) 2 

U a nn ' J - o T? 1 - 



app 'j g -gE /2a 



vv ' VV ' - exp (-gAE ,), or A , = exp(-gAE ,) 

11 r ^ , uv VV 11 ' VV 

V V V V vv 



Since A , is a probability, A , < 1 . 

vv vv 



<m> 



3S,nQ 



1 3Q 






3(eE a PP ) " « 



Evaluating the trace gives 



So if 



11 . -BE , 
v v ' 



e vv < 1 , we let A , = 



11 . -BAE , 
v v e vv and 



A , = 1 ; otherwise, we can switch indices, v <* v , i.e. 

v v 



* i -BAE , 



v v "i 
V V 



e v v and A ,= 1 . 

vv 



(b) To make a transition every move, we 



let A , = A , = 1 . So then 

vv v v 



w , = tt , 
vv vv 



p , V = v+1 
1-p, V = V“1 



! 

: : M ' ■ j 






4 

- i 

1 



/ dE 



3 /2 ° 2 cosh[g(A 2 + y 2 (E + E ) 2 ] 1/2 ] 

app ; 



2(A/ji) / dE e BE /2L cosh(gg) , and so 



3(BE Q ) ‘ 2 / dE e 8E sinh(gE) | 1 2(E + E )(p/A) 

app -co P C arm" ^ 



app 



= 2(A/y) / dE e~ 8E /2L 



sinh(g£) M (e + E app )/£ 



This form, which is independent of v, is possible because of the 
symmetry of the problem — each pair of adjacent levels is 



from which the desired formula follows directly. 

6.16 The modification to include windows is simple. Since we are basically 
adding a potential that is infinite outside the window (i.e., not under 
the umbrella) and zero otherwise, we need only check if we’ve stepped 



outside the window. This is done at the time we also check 




As described above, we determined T c for each value of L by 
extrapolating to M min = 0. The results of those extrapolations are 
given in the table. We then extrapolated T 0 to the L -*• » by using a 
procedure known as finite size scaling . In particular, we assumed 

[T c - T c (L)] b - 1/L 

where T o = T c (») . We determined the exponent b from the slope of 
C T o " T 0 (L) [! vs - 1111 L • This procedure gave b « 1.0 ± 0 . Then a plot 
of T c (L) vs. 1 /L gave 

T q = 2.26 ± 0.02 

as the 1 /L -> 0 intercept. This estimate can be compared with Onsager's 
T q = 2.269 . 



L 


T C (L) 


10 


2.60 


20 


2.43 


30 


2.38 


40 


2.34 



A discussion of the finite size scaling method is given in Sec. 1.2 of 
Binder and Stauffer's Chapter in Applications of the Monte Carlo Method 
in Statistical Physics , ed . K. Binder (Spr inger-Verlag, N.Y., 1984). 
Also see Part 2 of Gould and Tobochnik's Computer Simulation Methods 
referred to at the beginning of this manual chapter. 




Chapter 7. Classical Fluids 



7.10 The kinetic energy and pair potential parts of <E> will be the same as 
discussed in the text. For the three-body term, we write 

< l u (3) (r.- r r.- r) > - i N(N-1)(N-2) <u (3) (r r )> 
i > j > 5L= 1 _1 ~ J ' J ~ l 6 ~ 12 ~ 23 



N(N- 1 ) (N-2) 



/ dr N u (3) (r 12 ,r 23 ) e 



„N e "6U (r ) 



1 J dr 1 / dr 2 / dr 3 u (3) (r 12 , r^) N(N-1)(N-2) / dr N " 3 e" BU(r ] 



-0U(r N ) 



■ \ / ^ / di : 2 / d c 3 p (3/N) c 2 > c 3 > u(3)( Hi 2 > c 23 > 






For a uniform (homogeneous) system, p (3/N)( r ^, r ) 

3 ( 3 ) 

= p g (r 1 0 , r 2 3 ) , and we can integrate out one of the degrees of 
freedom to get the volume: 

= i V / dr 12 / dr 23 p 3 g (3) (r 12 , r^J u (3) (r_ 12 , r^) 

“ N • 5 ^ / ^-1 2 / ^23 g ( 3 )( r - 12 ’ ^ 23 } u ( 3 )( - 12 ’ C 23 } ’ 



= | k B T + \ P / d 1 8(h)u(r) 



6 p / d -i 2 ^ ^23 8 ( -12’ -23 ) U ( -12’ -23 ) 



7 ,, Qn fS&nCh 3 r N -0U(r ) 

7.H BP = l— i XiN = w in j dr e 



_ 1 / O 

We change coordinates to x. = V D r. , thereby accounting for each of 




the three cartesian components, dx^ = Vdr^ . Hence, 



Bp - an V N / dx N exp[-g J u(V 1/3 x )] 

i . n J 1 



where the integration is now over the unit volume, not the system 
volume. Differentiation is now easy, directly affecting the pair 
potential. Accounting for the sum over N(N-1)/2 equivalent pairs, 
differentiation yields 



J dx N N(N-1)exp[-B l u(V 1/3 x )][du(V 1/3 x 10 )/d(V 1/3 x )] x V 1/3 
i>j=1 12 12 J 1 2 V 

/ dx N exp[-g l u(V 1/3 x. .)] 
i> j = 1 1J 



P - y /dr N N(N-1)exp[- B J u(r..)][du(r 12 )/dr 12 ]r 12 /J dr N exp[- B f u(r .)] 



h I* 1 J a C 2 r 12 »(»-:) / •Sr' 1 ' 2 e"*“ (r ") / j ar » e -»U(r") 



’-hi*, r ,2- <2/N> <c,.c2) 



Therefore, for a uniform system 



1 - (Bp/6) j dr g(r)r du(r)/dr . 



12 / dr[e- Su(r) - i] . f dr r 2 [e -Bu(r)_ ^ 






As r ■* » t u(r) ■* 0 , and 

rl''"-']. ^ 0 



for quickly enough decaying u(r), i.e., faster than r“3. Therefore, the 



boundary term in the integration by parts gives zero. As a result, 



7.13 



7.1^ 



i r -Mr) 1 f J 



f dr[e 



- 1] - j i- 4) / dr % 
2 0 



2 -Bu(r) du 
r e r d? 



Bu(r) du 



- - 8/6 J dr e‘ BUUW r 



dr 



For hard spheres 



B 9 (T) = - 2 ir[- J dr r 2 ] = ~~ a b . 



0 



For a square well 



,(T) = - 2 ir { - / dr r 2 + / dr r 2 [e® e - 1 ] } 



2 ir r >3 „B£/ '3 3 n 

= — La -e (a -a)]. 



At the Boyle temperature, 



p [o' 3 - e B B e (o '3_ ,3)] . o , 



therefore , 



e/k„T„ - in 



<3 3 

o 0 



T b „„ r_'3/c_'3_ ,3 



k^n (a°/(0 ,J - 0 J )) 



Q.E.D. 



The figure shows a graph of the square well B 2 with 0 ' adjusted such 
that 



2 r 

J dr r u LJ (r) = - e J dr 



which implies o' = 9 1 /3 o . With this choice of 0 ' , the square well 




7.28 The velocity or momentum averages are isotropic and factorable from the 
coordinates. Therefore, only one average is necessary to answer all 
parts to this Exercise. 

(a) <v 2 > = j dp v 2 ~ 6p /2m / J dp e~ 6p /2ra 

“ 2 -6p 2 /2m 

J d P x ( P x /m ) e 

« 2 

f . -Bp /2m 
J dp x e OH x 

— 00 

since the p^ degree of freedom is statistically uncorrelated from 
the others. Performing the integral gives 

, 2 1 

<v > = — — . 

x Bm 

Alternatively, 

2 2 1 2 2 1 2 

<y> = - -r m <v > = - - <p /2m> 




7.29 



(b) 



(c) 



(d) 



(e) 



2 2 
<v v > 
x y 



<v 2 Xv 2 > = (d- ) 
x y 1 Bm ' 



2 2 
<v > = <v 

x 



<v > = 0 
x 



2 2. 

v + v > 

y z 



3<v 2 > = J- 
x gra 



<(v + bv ) 2 > = <v 2 > + 2b<v > < v > + b 2 <v 2 > 

x y x x y y 



(Hb 2 ) <v 2 > 



1 

JN 



« -Bk(x - x ) /2 
[( / dx e 0 



2 

“ -Bp /2m 

M / dp e x )] 3N 



„3N 



2tt 

BK 



2m 7T -,3N 

D -J 



= (BKw) 3N . 



(c) Since e ® 10/12 - 1 + gMto/2 for small B , 



QM 



Ed + ^ ) - (l - ^“)]‘ 



1 +b 
Bm 



Classically, these averages won't be affected by isothermal 
volume/ pressure changes (i.e. through the potential) since the 
momentum distribution is dependent only on the temperature. When 
the classical assumption 'breaks down at extremely high pressures, 
one would then need to consider the quantum dispersion of particles 
since the dispersion (i.e., the uncertainty principle) leads to 
statistical correlations between positions and momenta. 

(a) Since the Hamiltonian is separable, one can factor the partition 
function as discussed in Chapter 4. This gives 



-BWm/2 



1 - e 



-BMw 



j 3N _ ^ e BIW2 _ e -BHw/2i 



-3N 



where we have used the fact that + n) Mw with w = /k/m is the 
nth energy level for a one-dimensional harmonic oscillator. 



(b) In the classical world, the partition function still factors, 
Q = 



h 



1 / dr N / dp N e^’ P N) 



( Bhu) 



-3N 



classical 



I 

I 7-31 (a) For p ■> p cp (but p < p Cp ) we approach the g(x) for a close packed 

system, i.e., delta functions spaced l apart. The lower the 
density, p, the broader the peaks, and the more rapid the decay to 
unity at large x. 



g(x) at high p 




x/l 



3N 



72 




Since u(x) is » for |x| < 8 and zero otherwise, g(x) is therefore a 

_ Qy ( X ) 

step function at low p. As p -> 0 , g(x) = e 



9 i 



g(x) at low p 




Notice for a hard potential system, 0u(x) is either infinite or 
zero, so gu(x) is independent of temperature. Thus g(x) is 
independent of temperature. The configurational properties of hard 
core potential systems have no temperature. 



fg» 



i 



i 

j 



2_ 

Ttmg 



and 



(Ul) < 1 mv 2 > - ^ - lk B T 



Finally, 

(iv) <E> = <KE> + <PE> 

■ N ‘ < i > - i N V 



where we have noted that <PE> = 0 since the "potential energy" 
is zero for all acceptable configurations. 



(f ) Since Bp = p + bp + . . . , we can identify b as the second virial 
coefficient. Hence 

.8 



-0u(x) 



b = - ~ j dx[e 1] 

— CO 

Note that b is independent of temperature. 



/ dx(- 1 ) = l . 
- 8 . 



Since the integral of pg(x) determines the number of neighbors, and 
since there will be one neighbor between 0 and (3/2)8,, we have 



1 



( 3 / 2 ) 8 . 



dx g cp (x) = p^J 



= l 



Note that the momenta are uncoupled from the positions, and that 

2 

the momentum distribution is proportional to exp(-gp /2m) . Thus, 

(i ) <v> = 0 . 

Further , 



(ii) < | v | > =2 / 



-gp /2m 



dp / / 



-0p 2 /2m , 
e dp 



! 

1 

! 

■J 



(g) The partition function for the system can be calculated exactly. 

The easiest way is to think of the excluded volume due to the rods 
(NS.) as being subtracted out of the possible position space of each 
rod (L). Secondly, configurations which are identical except for 
the switching of rods are counted as the same state, i.e., we need 
to divide by a factor of N! . So, 



Q = 



\_ 1 _ 
, N N! 



L-N8, L-N8, L-NS, 

/ d Xl / dx 2 ... / dx^j 



if Jrr (L “ m)N 



This value for Q can also be arrived at from scaling and changes of 




coordinates of the configuration integral: First, adopting one 

particular ordering of the rods, 



Q “ Tj / dx i J dx 2 "• J dx N e " 

X 0 " " 



. L-(N-1)J, L-(N-2)S, L 

1 J dx ! / dx ... j 

X l 'x + i x N-1 + SI 



which becomes after a change in variables to x* - x^ - U , 
L-NH L-NJ, L-NJ, 

« = ^ / dx ; ! dx 2 ••• / , s 

A 0 X 1 X N-! 

1 1 ( t „ t iN 

" .N N! (L NS,) ' 

A 



„ 38.nQ 1 p 

8P ~ 3L T N N L-NSt “ 1-pA 



Thus, b = a , as already noted. 



7.32 (a) Neglecting the internal structure in the low density limit, the 

gaseous A r is an ideal structureless gas. Thus, for M atoms 



(id) 3 anQ ' 



„ w „ r V ,2™. 3/2-, 
= JlnM - Sin (— 5 - (— 5 —) J 

h 3 8 



with p = M/V 



* P ( gs ) 3/2 ♦ *np 



= f(B) + 8 -np 



(b) Let U (R ) stand for the total potential energy for all N water 
w 

molecules with R N denoting all the coordinates necessary to 



describe the configurations of all N water molecules. In the 
density limit, different argon molecules do not interact with 



other, but they do interact with the water molecules. Therefore, 
quite generally, we can view the total potential energy of the 



solution as 



U (R ) + I U (r.j R N ) 
1 = 1 ~ 1 



where r , is the position of the i th argon atom and there are 
atoms. The partition function is then 



M such 



exp[-Nf w (B)-Mf(g)] ij- -L- /dR N /dr, . . .dr M exp{-B[U w (R N ) + l U(r.;R N )]} 



n e" MF(8) , N -BU w (R N ) . M - 8 U(r.iR N ) -0U (R N ) 

Q w ~ M! / dR e C/ d r 1 .../<lr M H e ]/JdR N e w 

1 = 1 



-Mf(e) 



/dr 1 .../ 



M -BU(r. : R") 



dr.. < n e 
~ M i = l 



where Q w and <...> w denote the partition function and ensemble 
average for the pure liquid water solvent. Now, since the argon is 
dilute, different argon atoms are uncorrelated. Thus, the average 
of the product can be replaced with the product of the average. 
Hence, due to indistinguishability , 

-Mf ( B ) -BU(r ;R N ) M 

Q " V M <e _1 > 

w M! w 



_ n „M e' MCf<8) + 8A ' j] ba 11 _8U( Cl’ Rl,) 

- Q V — , where e 8Ay = <e 1 > . 



Note that since the solvent is isotropic, the last average is 



independent of r 
-i 



The desired formula for y is obtained on differentiation of 
JlnQ with respect to M with Q given by the last equation. 



Due to <1> w =i, it is obvious that Ay = 0 when U(r 1 ; R N ) is zero. 



(o) At phase equilibrium, the chemical potentials are equal. Hence, 



f(B) + inp G = f(B) + S,np L + gAy , implying p Q = p^e 



BAy 



dAy. 



P„ l dr ~ ‘W'~;A) u flU (r) . 



dX K w 1 ~ AW ’ “AW' 

Integration from X = 0 to X = 1 gives the final desired result. 



and since for an ideal gas, Bp = p , 

G 



Bp - P L e 



BAy 



Q.E.D. 



(d) We can simply plug into the result for Ay derived in Part b. 
Alternatively, we can follow the suggestion noting 



Ay = change in Helmholtz free energy due to introducing 
one argon atom into pure liquid water. 

= - k T S,n[Q( 1 argon, N waters )/Q ] 

D W 



7-33 To prove the virial theorem in two dimensions, use the same proof as 



r . On 



given for Exercise 7.11, except that now, ^ = V~ 1/2 
differentiation, this brings down a factor of 1/2 rather than 1/3. One 
therefore finds 



Bp/p - 1 - (Bp/M) / dr g(r)r ~ ^ r ^ 

= 1 - (Bpu/2) J dr r 2 g(r) du ^ 



dr 



from which the desired result for exp(-gAy) follows immediately. 

Next, define Ay according to 
A 



-BAy. 



e = < n e 

i 



Then 



~ BAU AW ( k - !l7l> 



. -BAu (|r-r |) 

d B Ay x <® MU AW C l=-^1 I > J « > w 

dX = " BXU AW ( l r " r 'il ) 

< n e AW 1 > 



/® N N u flW (k-Cil) 



-BU, 



JdR N 



-BU, 



3 B<N U AW>X 



where 



U X 3 u w (R 5 + .1 x U AW ( I ) 



and <...> x indicates the ensemble average weighted by that 
X-potential energy. Thus, from the definition of radial 
distribution functions, 






4 



■ 



1 + (pir/2) / dr r 2 [y(r) e Su ( r ’)] 
0 ar 



where we have introduced y(r) defined so that g(r) = exp[-Bu(r )]y (r ) . 
For hard disks, 



-Bu(r) 

e = 0 , r < a 

- 1 , r > a . 



Thus , 

d -gu(r ) 
dr e 

As a result, 



6(r-o) 



Bp/p = 1 + pn/2 / dr r 2 y(r) S(r-ct) 
0 

= 1 + (piro 2 /2) y(o) 



Recall that y(r) is the distribution function for cavities. It should 

be continuous at r = a . However for r > a, g(r) = y(r) so 

y(a) = lim g(r) . Hence, 
r+a + 

Bp/p = 1 + (po 2 tt/2) g(o ) 



7-35 



(a) Let U = l [u (r ) + Xu (r )] . Then from 

A U 1J I 1J 



dA/dX = - 



1 dilnQ 1 d 



0 f - “ BU x 

In j dr e 



B dX B dX 

one finds using arguments similar to those in Exercise 7.32 that 
dA/dX = (p 2 V/2) / dr g^CrJu^r) . 

Integration over X from 0 to 1 yields the desired result. 



I 



8.7 



f N -» ( V V r N " BU o 

(b) Q/Q 0 = Jdr N e / /dr e 



-eu n -b<u 1 > 0 

“ < e >0 * 6 



Further , 



<u i>o = ‘ < 



N(N-1 ) 



u, (r . ) > r 



= (p 2 V/2) / dr g Q (r)u 1 (r) . 



Therefore, since A - A„ - - k T An(Q/Q ) , 
U D U 



A-A 



0 < j p /dr g Q (r)u 1 (r) 



£ 



■ 



Chapter 8. Statistical Mechanics of Non-Equilibrium Systems 
Since H [z] =1 z > q* 

B i, 

0 z < q* 

and S[q(0) - q*] requires q(0) - q* , 
we interpret 

k BA (0) = k BA (0+) = k BA (e) ’ 

In other words, we're looking for the flux through q* for small enough 
time that trajectories which are started at q* essentially do not 
recross. Note the existence of v(0) guarantees an e small enough that 
q(e) = q(0) + ev(0) + 0(e 2 ) . Then 

k (0) - x 1 <v(0) < 5 [ Q ( 0 ) - q*] H [ q ( 0 ) + ev(0)] > 

BA A B £ ^ 0 + 

= x 1 <v(0) 6[q(0) - q*] H B [q* + ev(0)] > 

A e->0 + 

= x^ 1 <v(0) 6 C q ( 0 ) - q*] 0[v(O)] > 
where 0 is the Heaviside step function. We can then write 

k BA (0) - x" 1 <v(0) 9[v(0)] , 6[ _ q _ ( _ 0 ) - _q * ]> <fi[ (0) _ . 

BA A <S[q(0) - q* ]> 

* 

The first term is a conditional probability given q(0) = q . However, 
since q(0) is statistically independent of v(0) in the ensemble, we can 
drop the condition. Note that the original expression can't be factored 
in this way since v(t) = v[t ; q(0) , v (0) ] . 

So then 

k BA (0) = x' 1 <v(0)e[v(0)]> <6[q(0) - q*]> , 



and since the distribution for v(0) is even and the time origin is 



arbitrary, 



Since 



k (0) = < | v I > <6(q-q*)>/2x . 

dA ' 1 A 



( TST) 

Recognizing [q(t)] = 0[v(O)] gives the transition state theory 



expression for k 



(TST) 

B 



D = 1 / dt <v(0) -v(t)> 

i 0 



then 



3D % m 

t = *-p- = BmD 

<\ mv > 



1 d AD 2 A(-N 

6 dt AR (t) 



t+« 



1 ^ 2 S 
2 <V > T 



8.13 |r AR 2 (t) =2 / dt '<v(0) • v(t ' )> =2 / dt'<v 2 > e' t,/T 



= 2<v 2 > t [ 1 - e t/T ] . 



For Ar, MW = 40 g/mole, so at T = 80K 
t = 60 fsec . 



Integration and enforcing the initial condition of AR (0) = 0 gives 



AR 2 (t) = 2<v 2 >tt + 2<v 2 >t 2 [e t/x - 1 ] . 



For t << t , to 2nd order in t. 



AR 2 (t) = 2<v 2 >xt + 2<v 2 >t 2 [1 - - + ~ - 1] 

T 2t 2 



2 2 ? 

= <v >t £ = <(v Q t) > , 



i.e., inertial (non-diff usi ve ) behavior. For t » t , 



AR 2 (t) ■* 2<v 2 >tt - 2<v 2 >t 2 = 2<v 2 >x(t-T) 



which represents the long time diffusive behavior. The function is 
plotted for all time in the figure. 




■ 



} 



8.20 The primary variable x(t) drives the bath, and the driven bath affects 
the primary variable through the force 

CO 

f(t) = f (t) + / dt' x h (t-t') x(t') 



■ f b (t) - 6 f m dt ' c c b (t " t,)] x(fc,) 



= f b (t) + ec b (0) x(t) - 8C b (t) x(0) 



/ dt' C b (t-t’) x(t') (a) 

th 

C b (t-t’) , and the third arises from integration by parts and the fact 



where the second equality uses the connection between x b (t-t') and 



(essentially a notational convention) that f(0) =f=Ycv =f(0) 

. i i b 

l 

f b (t) differs from f(t) only as time evolves from the initial phase 
space point at t = 0. 

Now consider the average of f b (t), averaged with the initial 
conditions x and x held fixed: 



<f b (t)> 



x,x fixed 
at t=0 



< f b (t ; x , x , {y A ,y ± } ) 



e gxf > 



1 



b <e Bxf >. 



This equation follows from the fact that the distribution of initial 
bath variables differs by the factor (gxf) from that when the bath is 



82 



uncoupled from the primary variable. Let us employ the y^'s that are 
the normal modes of the bath, where mode y^ has the frequency u. . Then 



<f b (t)> 



x,x fixed i 



Bxc.y. 

I o <[y cosu.t + (y./iD. )sin<u.t] e 1 V — — L — 
jiii ill b Bxc.y 

<e 1 

v 3 . » XC i y i 



I c i 008 “ i t 38^7 * n <e 



= 8x l c 2 <y 2 > cos w i t = 8xC b (t) (b) 

where the second to last equality follows from the fact that each y. is 

an independent harmonic oscillator variable and therefore obeys a 

Gaussian distribution. The last equality is true since for normal mode" 
2 

<Y i Yj> - <Yj> <5.j , hence 



I c c <y y > cosU.t) = l c 2 <y 2 > cosio.t 
i,j J J 3 i 1 1 1 



By combining Eqs . (a) and (b) with mx(t) = f Q [x(t)] + f(t) and 
averaging over initial conditions of the bath variables, the generalize 



Langevin equation follows directly. 



: 



8.22 (a) <V(0)«V(t)> 



Argon Gas 



<y2> 



The very slow decay has a time scale of t , the mean collision ti 
in the gas. 





The oscillations persist for long times because for small 
displacements the atoms effectively sit in a harmonic well with 
period T . The slow decay is due to collisions with other atoms 
and also anharmonicities in the well which cause dephasing between 
different trajectories. 



<v^(0) v^(t)> 




The decay in the oscillations is due to the processes described 

above. Note that even if the lattice were one-dimensional, the 
2 2 

minima of <v (0)v (t)> would not drop, to zero because the ensemble 




average is over various initial conditions, so the turning point 
times are averaged over. In three dimensions, we have orbits with 




8.23 



no turning points. This is reflected in the three separate phases 
that are averaged over: 



<v " (0)v 2 (t )> = <(v 2 (0) + v 2 ( 0) + v 2 (0) ) ( v^{ t ) + v 2 (t) + v 2 (t))> 



= 3<v 2 (0) v 2 (t)> in the case of an 
isotropic lattice. 



(d) 



<u(0)- u{t)> 




The small oscillations are due to the vibrations in the lattice. 
The long time decay is due to the infrequent process of a CO 
molecule flipping over completely in the lattice. 



For a small perturbation from the equilibrium zero-field distribution, 
the fluctuation-dissipation theorem tells us the rate of relaxation to 
equilibrium is the same as that of a spontaneous fluctuation to that 
non-equilibrium distribution. At time t = 0, the field is shut off but 
the flow v(0) * 0 . Since the flow will dissipate with time, 



v(t) = Av(t) a <6v(0) • 6v (t )> = <v(0) *v(t )> 



v(t) 
v (0) 



<v(0) *v(t)> 
<v 2 > 



The relaxation time is then 
^ <v(0) -v(t )> 

0 



relax 



= J dt 



<v 2 > 



; 

: 

■f 

f 

1 



0r T relax = ® mD 



.26 Substituting c^(t) » c - c^Ct) - c^Ct) into the rate equations gives 



o^t) = - 0< 31 + k^) c^t) - k 13 c 2 (t) + k^c 



c 2 (t) = - k 23 ( t ) - (k^ 2 + k 23 ) c 2 (t) + k^c 



or, in matrix notation, 



d_ 

dt 



k 31 + k 1 3 
k 23 



13 

k 32 + k 23 



13 

k 23 



in terms of the 



The homogeneous part can be solved by expressing 
eigenvectors of the above matrix, while the inhomogeneous part just adds 
a constant to d 1 and d 2 : written in the basis that diagonalizes the 

above matrix, 



dt 



d 1 




X 1 


\ 

0 




1 \ 
d 1 


d 2 




0 


A 2 




d 2 



d 1 (t) = f ^ 



•Alt 



const 



-x 2 t 

d 2 (t) = f 2 e + const' 



where f^ and f 2 are constants determined by the initial conditions. 
Thus the concentrations c^ and c 2 can be written as the sum of two 
decaying exponentials, 

- A. t 



(t ) = A 1 e 



-A 2 t 

+ B.j e + < 0 ^ 



-At -At 

°2 (t ) = A^e + B^e + <c 2 > 



where „ and 



(a) So, 



are eigenvectors of the above matrix. 



-A.t -At 

Ac 1 (t ) <= e + B e 



-A.t -At 

Ao 2 (t) = A 2 e + B 2 e 



where 



A 1,2 = CCk 31 + k 13 5 + (k 32 + k 23 )]/2 



± /7T( V k 13 ) - (k 32+ k 23 >]/2 r + k 13 k 23 



(b) For e << 1 , then ^ 3 , k 23 » k , k 32 . 



We can write 



A 1,2 = [(k 31 + k 32 } + (k 13 + k 23 )]/2 



+ / (C0< 31 + k^) + (k 32 + k 23 ) ]/2 } (k 3i k 32 + k 31 k 23 + k 13 k 32 ) 



C(k 31 + k 32 )/2 + (k i 3 + k 23 )/2] 



x [1 ± A - (k 31 k 32 7 + k 13 k 32 )/[(k 31+ k 13 ♦ k 32+ k 23 )/2] 2 

A 1 2 is now conveniently written in terms of the small parameters 
k 31 ’ k 32 ’ so that we can ex P and the root /1-e = 1 - 2 e : 



A 1 , 2 " CCk 31 + k 32 )/2 + Ck 13 + k 23 )/2] 



X {1 ± [ 1 - 1 (k 31 k 32+ k 31 k 23 ♦ k 13 k 32 )/((k 31+ k 13 ♦ k 32+ k 33 )/2} 2 ]} 



So the + root is 



[( k 13 + k 23^ /2 + C2 - 0(e)] = k 13 + k 23 - T transient 



and the - root is 



(k 31 k 32 + k 31 k 23 + k 13 k 32 )/(k 31 + k 13 + k 32 + k 23 } 
“ ^ k 3i k 23 + k 13 k 32^ / ^ k l 3 + k 23 ) = T rxn ‘ 



From the above , 



-1 ( k 23 



- (i 



rxn 31 



) + k ( ii ) « t _1 

1 K 32 k, _ + k„' transi 



'13 23 



ent ’ 



and so the relaxation is dominated by T rxn . When = k i 3 



t = (k + 
rxn 31 32 



(c) As shown above, the faster transient decay occurs on a time scale 
of 



(k 1 3 + k 23 } 1 ~ k 23 lf k 23 - k l 3 



(d) The two decay rates are analogous to the two rates in the reactive 



flux description: t 



mol 



t. . << T 

transient rxn 



The connection can 



be made by imagining preparing the system at the transition state, 
i.e., in state 3. Then the decay into states 1 and 2, c^t) and 
c 2 (t), follow the two decay rates, one much faster than the 
other. But c^t) <x <6n A (0)6n A (t )> by the regression hypothesis, 
and the time derivative of <6n A (0)6n A (t)> is just the flux in the 
reactive flux picture. In particular, t ' is on the order of k 31 



or k 32 , and is the plateau value for the reaction rate. 



(e) This is similar to the transition state theory idea where e - ^® is 
the probability of getting to the transition state 3 , and D « 1 /n 
is the rate to cross the barrier (3) once there. As we showed 
previously 

k (TST) <* < | v | > < 6 (q-q*)> . 

8.28 T = 300 K n = 0.01 poise D flr = 1 x 10~ 5 cm 2 /sec 

m Ar = 40 amu m H 0 = amu 



(a) and (b): 

The velocity distribution is just the Maxwell-Boltzmann distribution and 
is independent of the potential. So 



. 2 x 2 1 2 

< V > = v 2 



1.87 x 10^ cm 2 /sec 2 



or 

2 1/2 ii 

^ v > = 4.3 x 10 cm/sec 



for Ar atoms in the vapour and in solution. 



(c) Assuming inertial motion in the vapour, the time to move 1 0A is 

t = 10A/<v 2 r > 1/2 = 2.3 psec . 

2 

(d) Since AR (t) = 6 Dt , the time to diffuse 10A in solution is t = 170 



psec. We are justified in assuming diffusive motion since 170 psec 
is long compared to the velocity relaxation time calculated in part 

(e). 





